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We develop a theory for electron-electron interaction-induced many-body effects in three- 
dimensional Weyl or Dirac semimetals, including interaction corrections to the polarizability, elec¬ 
tron self-energy, and vertex function, up to second order in the effective fine structure constant of 
the Dirac material. These results are used to derive the higher-order ultraviolet renormalization 
of the Fermi velocity, effective coupling, and quasiparticle residue, revealing that the corrections 
to the renormalization group flows of both the velocity and coupling counteract the leading-order 
tendencies of velocity enhancement and coupling suppression at low energies. This in turn leads 
to the emergence of a critical coupling above which the interaction strength grows with decreas¬ 
ing energy scale. In addition, we identify a range of coupling strengths below the critical point in 
which the Fermi velocity varies non-monotonically as the low-energy, non-interacting fixed point is 
approached. Furthermore, we find that while the higher-order correction to the flow of the coupling 
is generally small compared to the leading order, the corresponding correction to the velocity flow 
carries an additional factor of the Dirac cone flavor number (the multiplicity of electron species, e.g. 
ground-state valley degeneracy arising from the band structure) relative to the leading-order result. 

Thus, for materials with a larger multiplicity, the regime of velocity non-monotonicity is reached for 
modest values of the coupling strength. This is in stark contrast to an approach based on a large-IV 
expansion or the random phase approximation (RPA), where higher-order corrections are strongly 
suppressed for larger values of the Dirac cone multiplicity. This suggests that perturbation theory in 
the coupling constant (i.e., the loop expansion) and the RPA/large-A expansion are complementary 
in the sense that they are applicable in different parameter regimes of the theory. We show how our 
results for the ultraviolet renormalization of quasiparticle properties can be tested experimentally 
through measurements of quantities such as the optical conductivity or dielectric function (with 
carrier density or temperature acting as the scale being varied to induce the running coupling). Al¬ 
though experiments typically access the finite-density regime, we show that our zero-density results 
still capture clear many-body signatures that should be visible at higher temperatures even in real 
systems with disorder and finite doping. 


I. INTRODUCTION 

Three-dimensional (3D) Weyl-Dirac semimetals, mate¬ 
rials that possess one or more Dirac points at which two 
bands in the bulk touch with a linear dispersion, E± (fc) = 
zkhvp\k\, have been of great interest in recent years. 
These materials may be thought of as three-dimensional 
analogs of graphene or as solid state (and nonrelativis- 
tic, i.e., the velocity of light is taken to be infinite so the 
bare Coulomb interaction is unretarded although there 
is a static background lattice dielectric constant reduc¬ 
ing the overall value of the Coulomb coupling) incarna¬ 
tions of relativistic quantum electrodynamics since they 
are described by quasi-relativistic fermions with Coulomb 
interactions at low energies. Dirac semimetals were first 
theorized long ago by Herring^ and by Abrikosov and Be- 
neslavskii^. More recently, these materials have attracted 
considerable attention not only because they host a sys¬ 
tem of interacting massless Dirac quasiparticles, but also 
due to the fact that they lie at the intersection of many 
interesting topological phases of matter. This observa¬ 
tion has driven a concerted effort to identify theoreti¬ 
cally materials that possess symmetry-protected Dirac 
points in the bulh^ - —. Different topological states can 
be realized by breaking certain symmetries; for exam¬ 


ple, breaking inversion symmetry lifts the Dirac cone de¬ 
generacy, giving rise to Weyl semimetals ) 3 ' 4 ' 6 ' 10 ! 11 mate¬ 
rials which exhibit anomalous transport properties^ - — 
and topologically protected Fermi arc states on their 
surface^ ’ 20 ! 21 . Over the last two years, numerous ex¬ 
perimental works have confirmed that certain materials 
such as NaaBi and CdaAs 2 are Dirac semimetals by di¬ 
rectly observing Dirac cones in the bulk 22-27 and Fermi 
arcs on the surface— Many additional works observ¬ 
ing and elucidating transport properties^ - — and phase 
transitions^ - — in Dirac semimetals have appeared. Even 
more recently, the first Weyl semimetals were discovered 
experimentally^ - —, with evidence for Fermi arcs^ and 
the unusual transport signatures they entail^ - —. The 
theory developed in this work applies equally well to both 
Dirac and Weyl systems [i.e., chiral linearly dispersing 
gapless 3D systems with filled (empty) valence (conduc¬ 
tion) bands touching at singular Dirac points], and we 
refer to both as Dirac systems from now on in this paper. 
The number of such Dirac points (e.g., from possible val¬ 
ley degeneracy in the system arising from band-structure 
effects) or, equivalently, the number of Dirac cones de¬ 
fines the number of flavors or multiplicity (N) giving the 
ground-state degeneracy (excluding the spin degeneracy 
for a Dirac material or due to there being two separate 
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cones per flavor in a Weyl material) which would play an 
important role in the theory. 

A crucial feature of interacting Dirac materials is the 
field-theoretic nrany-body renormalization of the quasi¬ 
particle properties, i.e., the Fermi velocity, the quasi¬ 
particle residue, and the effective Coulomb interaction 
strength are no longer pinned at their noninteracting 
band structure values but acquire a scale dependence. 
The strength of the renormalization depends on the di¬ 
mensionless ratio of the Coulomb interaction to the ki¬ 
netic energy as given by the effective fine-structure con¬ 
stant: 


Attvf Hvfk : 

where e is the electron charge, vf is the Fermi velocity, 
and re is the effective lattice dielectric constant of the 
material (re > 1 in general for solid state materials). Be¬ 
cause of the linear band approximation, any correlation 
function is dominated by high-momentum excitations as 
is manifested by an ultraviolet divergent dependence of 
these correlators on a cutoff scale A. Through a renor¬ 
malization of the quasiparticle properties, this ultravio¬ 
let cutoff dependence can be removed at the expense of 
introducing renormalized parameters which now depend 
on the scale (such as momentum, energy, temperature, or 
density) at which they are measured. This scale depen¬ 
dence is dictated by renormalization group (RG) equa¬ 
tions. For example, the coupling ID obeys 


da 
din fi 


( 2 ) 


where /r denotes the renormalization scale. The coeffi¬ 
cient on the right-hand side is called the beta function 
and is uniquely determined by the divergence structure 
of the theory. If a(fi o) is measured at one scale /io, 
Eq. ID then predicts the value of a(fi') at a different 
scale //. The quasiparticle renormalization is an observ¬ 
able effect and has been measured most notably in the 
two-dimesional (2D) Dirac material graphene in a wide 
variety of experiments using energy, density (which deter¬ 
mines the Fermi energy), and momentum as the running 
scaled— . 

In principle, temperature, provided it is much higher 
than the Fermi energy, could also be a varying scale to 
study the RG flow, and may very well be the most suit¬ 
able scaling variable for 3D Dirac systems in terms of 
experimental investigations. We note that our approach 
to the RG flow of quasiparticle renormalization does not 
make any explicit use of the precise value of the ultravio¬ 
let cutoff A since we express renormalization at one scale 
simply in terms of that at another value of the scaling 
parameter. This completely eliminates choosing an arbi¬ 
trary value of the ultraviolet cutoff which is theoretically 
ill defined except as an ultraviolet cutoff, making it prob¬ 
lematic to assign its precise value. In particular, the basic 
tenet of the RG approach is that physical quantities vary 
with the energy or momentum scale (and equivalently, 


with carrier density or temperature), but do not depend 
on the ultraviolet cutoff scale which should not show up 
in observable quantities. We should, however, mention 
that, in solid state systems, unlike in relativistic quan¬ 
tum field theories [e.g. quantum electrodynamics (QED) 
or quantum chromodynamics (QCD)], there is indeed a 
true ultraviolet momentum cutoff given by the inverse 
lattice constant. However, this cutoff is by no means a 
sharply defined unique quantity for calculating renormal¬ 
ized quasi-particle properties (e.g., it could be multiplied 
by 27r or some other constant or one could simply use 
the bandwidth as an ultraviolet energy cutoff), and thus 
any explicit calculation of quasi-particle renormalization 
using such a lattice (or band) cutoff momentum (or en¬ 
ergy) is not quantitatively meaningful in systems, such 
as Dirac materials, manifesting ultraviolet divergences. 

In this work, we present a calculation of the renormal¬ 
ization group (RG) equations for coupling, Fermi veloc¬ 
ity, and quasiparticle residue up to second order in the 
interaction strength a. Previous theoretical works have 
investigated renormalization phenomena in 3D Weyl or 
Dirac systems using the random phase approximation 
(RPA) or a large-A^ expansion 44 — as well as first- 
order perturbation theor y 13 : 41 : 64 . One of the primary 
motivations for our current work is to explore the effects 
of the leading corrections to RPA and to identify the pa¬ 
rameter regime in which a perturbative approach is valid. 
Similar investigations played an important role in the 
case of 2D Dirac materials, in particular graphene, where 
it was demonstrated that perturbative RG results are in¬ 
valid for all but very strongly screening media, and that 
apparent agreement between experiment and first-order 
theory is often completely spuriou s 65 : 66 . These findings 
led to the development of a theory that quantitatively 
describes renormalization effects in graphene that applies 
to a broader range of experimental setups^. It is there¬ 
fore of direct experimental relevance and of considerable 
fundamental interest to perform a similar analysis in the 
context of 3D Dirac materials to investigate how physi¬ 
cal properties of 3D Dirac materials vary under RG flows 
with higher-order corrections included. 

Therefore, in this work, we consider the effect of a 
Coulomb interaction on a system with N “flavors” of 
massless Dirac fermions to two-loop order. At this stage, 
N is simply a mathematical quantity defining the number 
of electron flavors or species, but in reality N is typically 
the (often large) valley degeneracy in the relevant 3D 
Dirac system. After reviewing the previous one-loop re¬ 
sults, we proceed to calculate the polarization, electron 
self-energy, and vertex function to two loops. We use 
the term “loop expansion” to mean the usual order-by¬ 
order diagrammatic perturbative many-body expansion 
in the Coulomb interaction, which is precisely equivalent 
to an expansion in the effective fine-structure constant 
for linearly dispersing Dirac systems. At one-loop order, 
we already see logarithmic divergences in both the po¬ 
larization and self-energy as a function of momentum. 
Therefore, at two loops, we expect to see log-squared di- 
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vergences, and we confirm this with explicit analytical 
calculations. We then consider the renormalization of 
the theory, in which we not only determine the necessary 
counterterms to cancel out the divergences, but also per¬ 
form several useful self-consistency checks on our results. 
We then make use of these results to derive RG equa¬ 
tions for the parameters in our theory, namely, the Fermi 
velocity vf, the effective fine-structure constant a, and 
the quasiparticle residue. 

While the remainder of this paper gives a detailed ac¬ 
count of both the systematics and the calculations of 
the renormalization group equations, we summarize here 
our main findings, which are the renormalization group 
equations obtained up to the two-loop expansion calcu¬ 
lation for residue, Fermi velocity, and coupling given in 
Eqs. (I196D . (12001) . and (12051) . respectively: 


d In Zjf, 
din /i 

din vf 
din fj, 


da 

din n 
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with C « 1.33318 a constant that we determined numer¬ 
ically [cf. Eq. ([72]) ]. We find that, in contrast to the 
one-loop result that implies a always going to zero as 
we go to low-momentum scales, the two-loop result re¬ 
veals that there is a critical value above which a instead 
diverges, just as in the case of graphene^. However, 
this critical value is much larger than in graphene; in 
3D Dirac materials, it is a c = 14.1298 (1 + -^-), which 
should be compared to a c = 0.78 for graphene. This 
would seem to imply that perturbation theory is reli¬ 
able over a much wider range of values of a in the 3D 
case relative to the 2D case. In particular, with exper¬ 
imental values of a ~ 0.1-1 in typical 3D Dirac mate¬ 
rials (e.g., in Cd 3 As 2 , dielectric constants in the range 
k ~ 20-40 have been measure d 67 : 68 , while velocities are 
typically^— in the range 10 5 -10 6 m/s), we might expect 
perturbation theory to give quantitatively reliable results 
since the condition a -C a c applies here. However, we see 
from the equation for the velocity renormalization that 
the second-order correction has again the opposite sign 
relative to the first-order result, producing a second spe¬ 
cial value of a, a* = 3n/N , where the second-order term 
is equal in magnitude to the leading-order term. It is 
generally the case that a* < a c , and for a* < a < a c , 
we find the unusual situation in which both the veloc¬ 
ity and effective coupling decrease as the energy scale 
is reduced. This trend continues until a reaches a*, at 
which point the velocity reverses its trend and begins to 
grow as the non-interacting fixed point in the infrared 
is approached. Moreover, we see that a* decreases with 
increasing Dirac cone multiplicity N, so that for larger 


values of N, a > a* is achieved for modest values of the 
interaction strength. This observation may be particu¬ 
larly relevant for the pyrochlore iridates where N = 13^, 
meaning that for a ~ 0.1, the second-order term in the 
velocity constitutes a large 25% correction. On the other 
hand, for CdaAs 2 with N = 1, this term yields a few- 
percent correction to the leading order. If we interpret 
the appearance of a* as signifying a breakdown of pertur¬ 
bation theory (given that for this coupling, the second- 
order corrections are comparable to the leading order), 
then a* can be used as a criterion for determining the 
validity of results based on the perturbation theory in a. 
The fact that perturbation theory appears to work bet¬ 
ter for smaller degeneracies is particularly interesting in 
light of the fact that methods such as RPA or large -N ex¬ 
pansions work in exactly the opposite regime, suggesting 
that perturbation theory may be complementary to these 
approache s 62 : 66 . Thus, the perturbative loop expansion 
may be a reasonable approximation for 3D Dirac systems 
for small values of N, whereas in graphene it is only rea¬ 
sonable when the effective coupling constant itself is very 
small (for example, for substrates with very large dielec¬ 
tric constants so that effective k values are large). 

We also discuss how this predicted renormalization can 
be observed experimentally, namely, by measuring the 
plasmon frequency and the optical conductivity. While 
the control over the doping in experiments on a 3D 
Dirac material is more restricted than in experiments 
on graphene (since gating is not an option for 3D ma¬ 
terials as it is for 2D graphene), the temperature of the 
system can serve as another energy scale. If this scale is 
much greater than the Fermi energy, which is determined 
by the doping density, then we expect that our results, 
which have been derived for an intrinsic (undoped) sys¬ 
tem, are valid for the system realized in experiments. In 
fact, having a finite temperature higher than the Fermi 
energy associated with the unintentional doping density 
in the system is a convenient and practical way of ap¬ 
proaching the intrinsic Dirac point physic s 63 : 69 : 70 . We 
elaborate on this argument by an explicit calculation of 
the Drude weight which is independent of initial doping 
even at intermediate temperatures of T > £f/2-, thus 
probing the intrinsic limit. While simple dimensional 
analysis predicts a linear temperature dependence of the 
intrinsic Drude weight, this linear scaling is violated by 
the logarithmic renormalization of the charge and Fermi 
velocity, giving rise to a superlinear temperature depen¬ 
dence instead. To leading order in a, the strength of 
this superlinear scaling is set by a renormalization group 
invariant quantity known as the Landau pole A l, a di¬ 
mensionally transmuted scale which marks the point of 
divergence of the one-loop fine-structure constant at high 
energy. To probe the relative effects of doping, temper¬ 
ature, and renormalization, we provide a full calculation 
of the plasmon frequency (which is related to the Drude 
weight) to leading order in the RPA, which clearly shows 
the transition from a finite-density extrinsic regime at low 
temperature to a zero-density intrinsic regime at high 
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temperature with superlinear logarithmic scaling viola¬ 
tions. An experimental observation of this super-linear 
temperature dependence will yield direct evidence for our 
predicted ultra-violet renormalization corrections. 

We summarize how our ultraviolet renormalization 
results compare with graphene, the 2D analog of the 
3D Dirac materials studied here, in Table U and with 
(3+l)-dimensional [(3+l)D] QED in Table ITT1 to which 
3D Dirac materials provide a nonrelativistic condensed 
matter analog. Let us first discuss the comparison 
to graphene, whose many-body renormalization effects 
were considered in detail in our earlier work s 65 ' 66 . We 
note that, while the Fermi velocity renormalizes in both 
graphene and 3D Dirac materials, the charge only renor¬ 
malizes in 3D Dirac materials (but not in graphene). This 
is because, in graphene (and, in fact, in all solid state 
2D Dirac systems), the Coulomb interaction is given by 
a non-analytic term in the action, and thus we never 
obtain any terms in perturbation theory that renormal¬ 
ize the charge. Physically, this is connected to the fact 
that the Coulomb interaction in graphene still has the 
3D Coulomb form since graphene is not the solid state 
analog of (2+l)-dimensional QED and is actually a 2D 
physical system existing in the 3D world. Therefore, the 
only renormalization of the overall interaction strength, 
as given by a, comes from renormalization of the Fermi 
velocity in graphene. On the other hand, the renor¬ 
malization of a in 3D Dirac materials is due to that of 
both charge and Fermi velocity [as in (3+l)-dimensional 
QED where both charge and mass manifest ultraviolet 
renormalization]; the Fermi velocity may be seen in this 
case as an independent parameter of the theory from a. 
One very important consequence of this is that, unlike in 
graphene, in which the system becoming weakly interact¬ 
ing at low energy scales necessarily means that the Fermi 
velocity diverges logarithmically at the Dirac point, the 
Fermi velocity can (and, as it turns out, does) actually 
remain finite even at low energies in 3D Dirac systems, 
since the non-interacting limit can be achieved by the 
charge renormalizing to zero at low energy. This gives us 
a very major distinction between the 2D and 3D Dirac 
systems: in 2D, we must see a logarithmically divergent 
Fermi velocity at low energy because it is the only way for 
the system to become weakly interactin g 56 ' 60 ' 65 , but this 
is not necessary in 3D. Note also that we expect to see 
a logarithmic divergence appearing for the vertex renor¬ 
malization in graphene at 0(a 2 ), even though it has not 
been calculated. This is because, at second order, the 
self-energy acquires a temporal component^, just as it 
does in 3D Dirac materials. Due to gauge invariance, we 
expect that the vertex renormalization should diverge as 
well at orders higher than the first, and in such a way that 
the Coulomb field strength does not renormalize at any 
order, as we see happens in 3D Dirac materials. Thus, 
as summarized in Table [H there are major mathemati¬ 
cal and physical differences between many-body effects 
induced by the ultraviolet renormalization in 3D and 2D 
Dirac systems arising in condensed matter physics. 


We now consider the comparison between 3D Dirac 
materials and (3+l)D QED, as summarized in Table ITT1 
Here, the situation is in some sense opposite to that of 
graphene. The interaction strength in 3+1D QED renor¬ 
malizes only because charge renormalizes, since the speed 
of light cannot renormalize due to the requirement that it 
be constant. Both, however, still exhibit the well-known 
Landau pole (as does graphene^), a divergence of a at 
a large, but finite, energy scale at one-loop order. While 
the Fermi velocity renormalizes (and, in fact, diverges 
itself at the Landau pole), it does not diverge rapidly 
enough to cancel the effect of the diverging charge. We 
also note that logarithmic divergences start appearing 
at lower orders in a in QED than in 3D Dirac materi¬ 
als in the vertex correction. This is because we treat 
a non-relativistic Coulomb interaction here, rather than 
the full relativistic electromagnetic force, as is done in 
QED. This means that the quasiparticle residue, corre¬ 
sponding to a temporal component of the self-energy, al¬ 
ready differs from unity at 0(a) in QED; in the 3D Dirac 
case, this only happens at 0(a 2 ). In both cases, how¬ 
ever, the vertex function and the quasiparticle residue 
are renormalized in such a way that the strength of the 
electromagnetic field remains unchanged, as is required 
by gauge invariance. This divergence of the vertex in 
QED at 0(a) also, ironically, leads to the cancellation of 
a potential log 2 divergence at 0(a 2 ). We will find that, 
for the 3D Dirac material case considered here, only the 
diagram in which the self-energy is inserted into one of 
the electron lines in the bubble has a log 2 divergence. 
In QED, on the other hand, the vertex correction con¬ 
tribution also has a log 2 divergenc e 73 ' 74 ; this divergence 
exactly cancels that from the self-energy correction con¬ 
tribution, leaving only a simple logarithmic divergence at 
0(a 2 ). 

One fundamental difference between solid state Dirac 
materials and QED is of course the large difference 
in the applicable effective coupling in the two cases. 
Whereas the fine-structure constant for QED is always 
j- sa 1/137, the corresponding effective fine-structure 
constant in Dirac systems (either 2D or 3D) is multi¬ 
plied by a factor of c/vfk [see Eq. (JTJ)] , which is typically 
very large since the Fermi velocity (~ 10 6 m/s) is much 
less than the speed of light (w 3x 10 8 m/s), although 
the presence of the effective background dielectric con¬ 
stant reduces the effective coupling in Dirac systems by 
a factor of order 10 typically in most materials (except 
for graphene suspended in vacuum where k = 1, giving 
a very large graphene fine-structure constant of around 
« 2.2). The largeness of the effective fine structure 
constant makes Dirac systems more like a strong-coupling 
QED problem (albeit with nonrelativistic Coulomb inter¬ 
actions) rather than the regular QED where the loop ex¬ 
pansion in powers of the fine structure constant is asymp¬ 
totic up to a very high (3> 137) order in the perturba¬ 
tion theory. We discuss the asymptotic nature of the 
loop expansion for 3D Dirac systems in this paper [see, 
e.g., Eq. (I204D ] and this issue was discussed in depth for 





5 



Grapheno 65 ' 71 

3D Dirac 

Vertex 

log divergence expected at O(o 2 ) 

log divergence at O(o 2 ) 

Self-energy 

log divergence at O(a), log 2 at 0(a J ) 

log divergence at O(a), log 2 at 0(a 2 ) 

Polarizability 

log divergence at 0(a 2 ) 

log divergence at O(a), log 2 at O(o 2 ) 

Charge 

No renormalization 

Renormalizes 

Fermi velocity 

Renormalizes 

Renormalizes 

Optical conductivity 

log divergence at 0(a 2 ) 

log divergence at O(a) 

TABLE I: Table summarizing the divergences appearing in graphene and 3D Dirac materials for various quantities. 


3D Dirac 

3+1D QEDS-Zi 

Vertex 

log divergence at O(o 2 ) 

log divergence at O(a) 

Self-energy 

log divergence at O(a), log 2 at 0(a 2 ) 

log divergence at 0(a), log 2 expected at O(o 2 ) 

Polarizability 

log divergence at O(a), log 2 at 0(a 2 ) 

log divergence at O(a) 

Charge 

Renormalizes 

Renormalizes 

Fermi velocity/speed of light 

Renormalizes 

No renormalization 

Electron mass 

N/A 

Renormalizes 


TABLE II: Table summarizing the divergences appearing in 3D Dirac materials and in 3+1D QED for various quantities. 


graphene in Ref. [651 where it was found that the loop ex¬ 
pansion may already fail at the first order for graphene 
suspended in vacuum. One may wonder, given the large 
coupling constant in Dirac materials, whether some as¬ 
pects of its interaction physics resemble QCD rather than 
QED. It turns out that one specific aspect of interaction 
effects in Dirac materials (for both 2D and 3D) does in¬ 
deed have some superficial similarities with QCD. This is 
the behavior of the Dirac system when the effective cou¬ 
pling is larger than the critical coupling (a > a c ) so that 
the running coupling increases with decreasing energy, 
eventually diverging at the Dirac point. This appears 
similar to the QCD strong-coupling behavior (provided 
one starts with a sufficiently large coupling a > a c in 
the beginning for the Dirac system, whereas in QCD of 
course this behavior is generic because of the structure 
of the beta function itself). But, this similarity is some¬ 
what misleading because the beta functions in the two 
cases are fundamentally different, and the divergent run¬ 
ning coupling low-energy behavior of the Dirac system 
we find (for a > a c ) may very well be just an artifact 
arising from the failure of the perturbation theory. In 
addition, the Dirac problem neither has confinement nor 
true asymptotic freedom (as is obvious from the exis¬ 
tence of the Landau pole in the theory), which are two 
hallmarks of QCD. 

Another qualitative difference between QED and Dirac 
systems worth pointing out is the role of the Fermi en¬ 
ergy (i.e., a finite chemical potential away from the Dirac 
point) as well as temperature in Dirac materials, which 
provide additional physical parameters (doping density 
and temperature) for the experimental implementation 
of the RG flow in solid state materials. Either the finite 
Fermi energy (arising from any finite doping of the sys¬ 
tem) or temperature, by itself, can act as a low-energy 
cutoff for the RG flow just as the band width or the 
inverse lattice constant acts as a high-energy or high- 
momentum cutoff in solids. But, in the presence of a 


finite Fermi energy (temperature), the other parameter, 
i.e., temperature (Fermi energy), can be used as a prac¬ 
tical parameter to induce the RG flow as long as tem¬ 
perature (Fermi energy) is much greater than the Fermi 
energy (temperature). This can be very useful in ex¬ 
perimental measurements where the pure undoped Dirac 
limit (or zero temperature) can never be reached. Since 
the Fermi energy (i.e., the chemical potential) depends 
monotonically on the doping density, one can change the 
doping density (at fixed low temperature) or the temper¬ 
ature (at fixed low doping density) as the scaling param¬ 
eter to study the ultraviolet renormalization of Dirac ma¬ 
terials. By contrast, in purely relativistic field-theoretic 
problems (e.g., QED, QCD), the concepts of doping, fi¬ 
nite chemical potential (or Fermi energy), temperature, 
etc., simply do not apply, and the system is always an 
intrinsic system (in our sense) and can never be detuned 
from the precise critical point. 

It is also instructive to compare the many-body renor¬ 
malization properties of 3D Dirac systems with normal 
3D metals (or doped semiconductors) which are often the 
textbook systems for studying many-body effects aris¬ 
ing from long-range Coulomb interactions in solids^r— . 
Three-dimensional metals are defined by the sharp exis¬ 
tence of a finite Fermi surface and the associated Fermi 
liquid properties, with ultraviolet renormalization or di¬ 
vergences not playing any role whatsoever since the en¬ 
ergy dispersion in these cases is parabolic, and there is 
no zero-density critical point because of the existence of 
a large band gap separating the filled valence band from 
the partially filled conduction band with the finite Fermi 
energy. Thus, the key physics of interest in Dirac ma¬ 
terials, namely, the many-body interaction corrections 
associated with ultraviolet renormalization, is not rel¬ 
evant for ordinary metals at all. In fact, there exist 
detailed RG analyse s 79 ’ 80 of the stability of the Fermi 
surface in 3D and 2D metals starting from high momen¬ 
tum and systematically approaching the Fermi surface, 
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with the conclusion that, except for the possible occur¬ 
rence of interaction-driven superconductivity at high or¬ 
bital momentum channels (the so-called Kohn-Luttinger 
superconductivity^) at exponentially low temperatures, 
the Fermi surface is stable and there is no ultraviolet di¬ 
vergence (in the QED sense) anywhere for 3D (or 2D) 
metals. Thus, the key physics of ultraviolet renormaliza¬ 
tion and scale-dependent RG flow with logarithmic run¬ 
ning of the coupling are features not arising in normal 
metals as they do in Dirac systems. Of course, there 
are strong Fermi liquid renormalization effects in nor¬ 
mal metals, as discussed extensively in many standard 
many-body theory textbooks^ - — , and these effects are 
non-perturbative since the dimensionless interaction cou¬ 
pling strength in 3D metals (which depends explicitly on 
the metallic electron density and is not a simple effec¬ 
tive fine-structure constant as in Dirac systems) is large 
(~ 5 — 6 typically), but these many-body corrections (of¬ 
ten quantitatively large) are not indicative of any fun¬ 
damental ultraviolet divergence in the system with no 
logarithmic cutoff-dependent renormalization corrections 
arising anywhere in the problem. In fact, the divergence 
one worries about in metals is the infrared divergence as¬ 
sociated with the long-range Coulomb interaction (and 
not an ultraviolet divergence associated with large mo¬ 
menta) which is “cured” by considering a loop expansion 
in the dynamically screened Coulomb interaction (rather 
than the loop expansion in the bare Coulomb interaction 
we consider for the Dirac system in the current work), 
which provides a controlled perturbation theory for met¬ 
als at high electron density (which is equivalent to weak 
metallic interaction coupling). At metallic densities, such 
an expansion is only qualitatively valid (since there is no 
quantum phase transition in the problem), but a com¬ 
parison with experiment (or with quantum Monte Carlo 
simulations) shows that metallic many-body corrections 
calculated using an expansion in the screened Coulomb 
interaction works reasonably well even at metallic densi¬ 
ties where the system is far from being weakly coupled for 
reasons not very clear at this stage.— Density-dependent 
(but not cutoff-dependent) many-body corrections sim¬ 
ilar to the metallic quasiparticle renormalization effects 
also show up in doped extrinsic Dirac materials as non- 
ultraviolet non-logarithmic subleading corrections, but 
we ignore them completely in our current work since they 
are completely negligible compared with the ultraviolet 
logarithmic corrections as the Dirac point is approached. 

The rest of this paper is organized as follows: In Sec. 
m we describe the (Euclidean-time) action for our sys¬ 
tem and the associated Feynman rules. Section mu is 
dedicated to reviewing the one-loop results and deriving 
the two-loop polarization, while Sec. El provides a de¬ 
tailed calculation of the two-loop electron self-energy. In 
Sec. 0 we formulate the renormalized perturbation the¬ 
ory, determine the counterterms that arise therein, and 
derive the RG equations for our model. In Sec. IVI1 we 
show how our results can be tested in finite-temperature 
experiments on finite-density samples even though our 


results are derived at zero density. We present our con¬ 
clusions in Sec. I VIII We provide derivations of the one- 
loop polarization and electron self-energy, the two-loop 
vertex corrections, and the optical conductivity within 
the Drude-Boltzmann approximation in the Appendix, 
along with an alternate derivation of the RG equations 
using the Callan-Symanzik equation. The latter provides 
yet another consistency check on our results. 


II. MODEL 

The model that we will be using is that of a 3D system 
with N “flavors” (e.g., valleys) of Dirac electrons inter¬ 
acting via a Coulomb interaction. We note that the total 
degeneracy of the 3D system we consider is 2N including 
both spin and valleys. We will work with the following 
Euclidean-time action (setting h = 1): 

N r 

s = -53 dtd 3 R(lp a 'y 0 do'lpa+V F 1p a "/ l di'l/ja 
a=1 ’ 

+ <M,7Va) + 2^2 J dt d 3 R , (6) 

where the 7 matrices form a Clifford algebra, i.e., 
(y'b Y} = 2<W' for /j and v = 0,1,2, and 3 and 
ip = There is an implied sum on i in the above 

action from i = 1 to 3. The fields, ip, are four-component 
Grassmann spinors, and N is the number of valleys. In 
all, this model results in 2 N Dirac cones. The factor of 
2 comes from spin in a regular Dirac semimetal, or from 
the presence of a pair of Weyl nodes for each of the N 
flavors in a Weyl semimetal. 

Here and throughout this work, we will employ a 
quasi-relativistic notation, i.e., we define k = (ko,k), 
d A k = dko d 3 k, ft = fco7° + v F k ■ 7 , and k 2 = k^ + Vp\k\ 2 . 
The Feynman rules associated with this action are the 
following: 


Go(k) = j, (7) 



Here, the straight lines are electron propagators, while 
the wavy lines are the Coulomb propagators. 


III. POLARIZATION 

We now turn our attention to the calculation of the 
polarization. We will first quote the known first-order 
result, and then embark on a detailed calculation of the 
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FIG. 1: Leading-order polarization bubble, Ils(q). 


second-order result. A general feature of these contribu¬ 
tions at all orders is that they are ultraviolet divergent, 
and thus we introduce a momentum cutoff A into our 
theory. 




FIG. 2: Self-energy (left) and vertex (right) corrections to the 
polarization, IIs_B(g) and II\/(q), respectively. 


A. First order 


FIG. 3: One-loop electron self-energy, Si (q). 



The first-order, noninteracting, contribution to the po¬ 
larization, shown in Fig. [TJ is well knowns, so we simply 
quote the result here for the sake of completeness and 
leave the derivation to Appendix [A] If one computes the 
polarization using the given action in Eq. ©, then this 
contribution is 


Hs (q) 


kf _ Arin /k 2 + (2A+l) 2 ][^ + (2A-l) 2 ] 


24:Tt 2 vf 




(z 2 + l ) 2 
2A + y 

z 2 + (2A + y ) 2 1 


( 10 ) 


where z = v ^ and A = |4j-. The integral in this ex¬ 
pression can be done analytically, but the result is rather 
complicated and not particularly illuminating. We see 
that ns(g) is logarithmically divergent as A —► oo, the 
entire divergence coming from the first term. The loga¬ 
rithmic and finite terms, the latter of which we will need 
later, in ns(g) are 


n / v |g| 2a 

n B\q) ~ — 2 ~x—A/ 

9 3tt 


*(— 

\q\) 2 V 4 


• (ii) 


Note that this expression has an additional factor of 27V 
compared to that found in Ref. This is because said 
reference considers a single Dirac cone, while we consider 
a system with 27V flavors of Dirac cones. 


B. Second order 


We now consider the second-order correction. Two di¬ 
agrams, both depicted in Fig.[2j contribute to this correc¬ 
tion. Note that the first diagram contains the first-order 
electron self-energy Ei(g), shown in Fig. [3] as a subdia¬ 
gram. We will simply quote the result here, leaving the 
derivation to App. iBl 


Si (q) 




vfQ • 7 


( 12 ) 


1. “Self-energy” correction 


We will first determine the “self-energy” term, which 
we denote by IIg£;(g). The integral for this correction is 

n SE {q) = 2 N J ^ATV[ 7 °Go(A : + g )7 0 Go(fc)E 1 (fc)Go(fc)], 

(13) 

which includes a factor of 2 due to the symmetry of the 
diagram. We now substitute in the bare Green’s func¬ 
tions and self-energy and evaluate the trace. The trace 
that we must perform is 

t^VtVtV) 

= 4(2^° - 1 )(5^S ip - S pi 6 vp + 6 pp 6 vi ). (14) 


Using this formula, we may rewrite the integral as 


3 (h - h + h) 


TlMq) = l ^ N A 

3 Ulv 

where U, .Z 2 , and I 3 are integrals over fcg, given by 


(15) 


h = 2 v F [ ^ 


k 0 {k 0 + qo)\kf 


-00 2tt (fcg 4 - v 2 f |fc|2)2[(fc 0 T go) 2 + v 2 p\k + q\ 2 } 
\kU-v 2 F (\k\ + \k + m 


2 [q$ + v U\k\ + \ k + q\) 2 ] 2 


(16) 


h = 2 vi 


k ■ (k + q)\k\‘ 


r dko 

l-oo 27T (fc2 + v 2 f |fc| 2 ) 2 [(fc 0 + go) 2 + v 2 p\k + q\ 2 } 

1 


2 \k\[q 2 +v 2 (\k\ + \k + q\) 2 } 

1 

'\k+m+^ F (\k\ + \k+m 

_9o_ 

\k + q][q 2 + v U\k\ + \k + m 2 


k-(k + q), (17) 


I3 = VF 


k ■ (k + q) 


-00 2tt (fcg + v 2 F |fc|2)p 0 + <? 0 ) 2 + v 2 F \k + q ]' 1 
























|fc| + \k + q\ 


■k ■ (k + q). (18) 


2 \k\\k + q\lqZ + vU\k\ + \k + q\r 

If we now substitute all of these into Ils^q), we obtain 

^\k+m 


nss(9) = --^ vf N J 


d 3 k \k\[ql~v 2 F ( 


1 - 


W lq 2 o+v 2 F (\k\ + \k + q \) 2 } 2 

k ■ (k + q) 


\k\\k+q\ 


| + in ( 4 
3 \\k\. 


(19) 


We now switch to a prolate spheroidal coordinate system, 
with the z axis being the long axis of the spheroids. The 
system is defined as follows: 

k±,x = -|<f| sinh/xsinz/cos0, (20) 

fc_L y = — \q\ sinh/xsinz/sin0, (21) 

fc|| = - |g| (cosh/i cos — 1), (22) 

where 0 < /i < oo, 0 < v < 7r, and 0 < 6 < 2it. The 
Jacobian of this transformation is 

J = - \q\ 3 sinh/xsinz/(cosh 2 y — cos 2 v). (23) 


Some useful identities are 


|fc| = -|<f|(cosh/i — cos v), 
+ = -|g|(cosh/i + cosi/), 

and 


k ■ (k + q) = -|g| 2 (cosh 2 y + cos 2 v — 2). 


(26) 


Using these identities, we may also write the Jacobian 
as J = ||g| sinh/j,sini/|fc||fc + q\. Making this coordinate 
change, IIs^;^) becomes, after performing the (trivial) 
integral over 9 , 


ns_E(<?) =-2W° 2 N 

■n z g 


(cosh y — cos v) (z 2 — cosh" y) 

(. z 2 + cosh 2 y) 2 
2 1/ 2A 

X sin" v —|— in - 

9 6 \ cosh /i — cos v 


I dy / dv sinh/xsinz/ 
o Jo 
2 


(27) 


We now make another substitution, namely x = cosh y 
and y = cos v. Doing this, we obtain 


n SE (q) = f 

F 9 J i 

(Z 2 + X 2 ) 2 


dx 

2 

9 


dy 


In 


2A 

x-y 


.(28) 


We note now that this integral is divergent, and thus we 
must impose a momentum cutoff \k\ < A. In terms of x 
and y, this translates to imposing a cutoff of 2A + y on 
the x integral: 

2 lol 2 C 1 r 2 ^+v 
n SE{q) = — ~^—^a 2 N dy dx 

7T 2 9 J -1 J i 

(x - y)(z 2 - x 2 ) 


(z 2 + x 2 ) 2 

x (i -y 2 ) 


2A 


2 1 , 

7T + ^ l n 

9 d \x — y J 


(29) 


We may simplify this expression further by interchanging 
the order of the x and y integrals, 


r»l r*2\+y 

j y L 



n 2A-1 

dx —»■ 

/ dx 


J 1 


r2\+l 

+ 

/ dx 


J 2A—1 


dy 


>-1 
/*! 


dy , 


c-2A 


(30) 


and evaluating the y integral. Doing these, we find that 

dx ~r~Ti --7T7 ^t/i(x) 


(24) 

n sE(q) = — -^\-^— a 2 N 

7T 9 


/•2A+1 ,5 

(25) 

+ / dx 

J 2A-1 \ Z 


o 2A-1 


(z 2 + x 2 ) 2 


rh(x) 


(31) 


where 


fi(x) = 


1 


-(x 4 -6x 2 -3)lnf — 
2 v \ x — 1 


108 _ 

+ x(49 — 3x 2 ) + 24x In 


2A 


Yx 2 — l) _ 


(32) 


1 2{x) = sk {" 57 (w) +160 (w) 1 

- 132^^j (x 2 — 1) + (x — 1) 3 (29 x + 75) 


+ 12(x — l) 3 (x + 3) In 


2A 

x — 1 


(33) 


We may now extract the leading divergence of Ils^g). 
To do this, we take the derivative of the integrals above 
with respect to A, evaluate the result, and then expand 
in powers of A. A In 2 term would correspond to with 
a coefficient twice that of the In 2 term, while a simple In 
term would correspond to j. The first integral is the sole 
contributor to the divergence of ns^g); we find that it 
possesses both a In 2 and a In term. We thus find that 
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Z 

FIG. 4: Plot of the function f(z) appearing in Eq. (1341) . giving 
the dependence of the finite term in IIsB(g) as a function of 
z = qo/vF\q\- 


n SE{q) is given by 

nss(g) ~ 


\q) 2 2a 2 


9' 

+ 3 In 


27tt 2 


N 


kl 

2 {z 2 + 4) 


+ 1 


In 


2 kT 


- ^a 2 Nf(z) + Ol- 
ir 2 g 2 \X 


(34) 


The finite term, f{z ), could, in principle, be found ana¬ 
lytically, but the expression giving it is extremely compli¬ 
cated, and thus we determine it numerically. This term 
comes entirely from the integral involving ( x ), as that 
involving / 2 (a:) goes to zero as A —> oo. We present a 
plot of f{z) in Fig. [I] 


2. Vertex correction to polarization 


We now turn our attention to the vertex contribution, 
Ily(q). This contribution is given by 


/ ^4 1 ^ n 

J 


d 4 p 

( 2 W 


x Tr 




K + i o t + t o t 


k 2 (k + q) 2 (p + q) 


2 1 „2 


| k — p\ 2 
(35) 


+4S nV0 S pV0^ § pa _ 5 w 5 vo + ( 36 ) 


Using this, we find that the integral splits into nine terms: 


Ji = 


x 

J 2 = 

X 

Js = 

X 

J 4 = 

X 

J 5 = 

X 

J 6 = 

X 

J 7 = 

X 

J s = 

X 

J9 = 

X 


d 4 k 


d 4 p 


-VvavrNj ^ j ^ 

koPo(Po + qo){k 0 + q 0 ) 
k 2 p 2 {p + q) 2 {k + q) 2 ’ 

16 ™’ pW / A±J 

kp(Po + qo)p ■ jk + q) 
k 2 p 2 {p + q) 2 {k + q ) 2 


1 


d 4 p 
( 2 tt ) 4 |fc — p | 2 


\6iravpN J 

ko{k 0 + q 0 )i 

r d 4 k r d 4 p 1 

(2tt) 4 J (27r) 4 |^_p|2 

o-{p + q) 

k 2 p 2 {p + q) 

IGnaVpN J 

koPo{p + q 0 

2 {k + q ) 2 ’ 

r d 4 k f d 4 p 1 

(2tt) 4 J (2tt) 4 \k — p\ 2 

■{k + q) 

k 2 p 2 {p + q)‘ 

I 67 ravpN J 

Po{Po + qoV 

2 {k + q ) 2 ’ 

r d 4 k r d 4 P 1 

(27r) 4 J (2t r) 4 |jfe_£|2 

fc • (fc + q) 

k 2 p 2 {p + q) 
107 ravpN J 
{k 0 + 9o)(Pc 

2 (fc + g) 2 ’ 

r d 4 k j • d 4 p 1 

(2tt) 4 J (27 r ) 4 |fc —p| 2 
i + Qo)k-p 

* :V(p+q) 2 (k + tY ’ 

f ft [ ,f p u - 1 

J (2tt) 4 J (2.7 r) 4 |fc_ p| 2 

k ■ (k + q)p ■ {p + q) 

k 2 p 2 {p + q)' 
I 67 ravpN J 
k-(p + q){ 

2 {k + q) 2 ' 

r d 4 k r d 4 P 1 

(2tt) 4 J (2t r) 4 |jfe_£|2 

k + q) ■ p 


k 2 p 2 {p + q) 2 {k + q) 2 , 


— 16naVpN 


d 4 k 

( 2 n ) 4 


k-p(k + q)-{p + q) 


k 2 p 2 {p + q) 2 {k + q) 2 


d 4 p 1 
(2tt) 4 \k — p\ 2 


(37) 


(38) 


(39) 


(40) 


(41) 


(42) 


(43) 


(44) 


(45) 


Evaluating the frequency integrals is a straightforward, 
if tedious, exercise; once we have done so and simplified 
the result, we obtain 


We now evaluate the trace and perform the integrals on 
ko and po. The trace that we need to evaluate is 

Tr(7 0 7 At 7°7 I '7 0 7 P 7°7 <T ) 

= 45^° S p0 {2S u0 S a0 - 6 ™) - 4J p0 <^ 0 (<r ( V' T + 5°°5 vp ) 

_ 4 § pV 0 S p 0 ^ S cr 0 + § P<t 5 vO) 


n v (q) = 


X 

X 


— 47T aV F N 


d 3 k 

( 2 ^ 


d 3 p _1_ 

(27r) 3 |fc_p ]2 


1 

ko 2 + v 2 A\k\ + \k + q\) 2 H + v 2 F m + I p + q\) 2 ] 
{-qlQi+v 2 F Q2), (46) 
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where Q\ and Q 2 are given by 


Q 1 — 


Qi = 


+ 


k-p _ (■k + q)-p k ■ (p + q) {k + q) ■ (p + q) 
\k\\p\ \k + q[\p\ \k\\p + q\ \k + q\\p + q\ 

(47) 


(|fc| + \ k + q\)(\p\ + \p + q\) 

\k\\k + q\\p\\p + q\ 


[k ■ p\q\ 2 ~ {k ■ q)(p ■ q) 


(|fc | 2 + k ■ q — \k\\k + q\)(\p\ 2 + p-q- \p\\p + q\)}. 


(48) 


as follows. II v,e(q) is the “easy” part, which we can 
do analytically, while Ily^g) is the “hard” part, which 
cannot be determined completely analytically. 

The “easy” part comes from the ^ in Eq. flSH) . and 
is given by, after changing variables to x = coshp, and 
y = cos is, 


n v,E(q) 


M 2 « 2 Ar(T 2 2 t 2 \ 


(55) 


where 


To perform this integral, we introduce two sets of prolate 
spheroidal coordinates, one for k (p, is, 0) and one for p 
(p', is', O'). The quantities Q\ and Q 2 become 


Q 1 


Q 2 


4 

(cosh 2 p — cos 2 zz)(cosh 2 p' — cos 2 v') 
x [cosh p sin 2 v cosh p! sin 2 v' 

+ sinh p sin v cos v sinh p' sin v' cos v' cos(0 — 0')], (49) 
4|g | 2 cosh p cosh p' 

(cosh 2 p — cos 2 I'Xcosh 2 p' — cos 2 is') 
x [sin 2 ^sin 2 is' + sinh p sin is sinh p' sin is' cos(9 — 9’)]. 


(50) 


Furthermore, 


\k — p | 2 = -|g| 2 [cosh 2 p + cos 2 v + cosh 2 p' + cos 2 is' 

—2 — 2 cosh p cos is cosh p! cos v’ 

—2 sinh p sin is sinh p' sin is' cos (0 — 9')] 

= I f(p,v,6-,p' ,v' ,0'). (51) 


We thus obtain 


n v(q) = 


X 

X 


kl 


2 a 2 


g 2 256tt 4 

f 00 dp' I 


N 


dp 


dis 


dO 


dis' 

1 


t 27r , sinh p sin v sinh p' sin is' 
d9 


(.z 2 + cosh 2 p)(z 2 + cosh 2 p') 


f{p , ^,0;p',is',9') 

z 2 Q'\ - 


2 r>' Q 2 

\d \ 2 


(52) 


where Q\ = (cosh 2 p — cos 2 ^)(cosh 2 p! — cos 2 v')Qi. 

Note that the integrand only depends on 9 and 9' 
through their difference, 9 — 9' , and that it is a peri¬ 
odic function of both. As a result, we may “shift away” 
one of these variables, say, O' , thus making the integral 
over that variable trivial. The other integral, in this case 
over 9, can then be done with the aid of the formulas, 



(56) 

(57) 


Note that I\ is logarithmically divergent; we must there¬ 
fore impose a cutoff, as before. These integrals may be 
done exactly; the results are 


11 = In 
- 2 + 

1 2 = A In 


2 2 + (2A- l ) 2 


z 2 + 1 


tan 


+ (A+ 2) 114 


z 2 + ( 2A + 1 ) 2 


2A+ 1 


— tan 


2 2 + (2A — l ) 2 
2A — 1 


r z 2 + (2A + l) 2 ] _ 1 z 2 + 1 — 4A 2 


z 2 + (2A- 1) 2 _ 


x 


tan 


-1 


2A + 1 


— tan 


2 z 
2 A- 1 


(58) 


(59) 


The divergence of II v,e(q) comes from I\. To determine 
the full divergence, we need both the In and finite terms 
of Ip. 


Ji “ 21 "(i) +1 "(tt?)- (601 

We can also determine the finite term; in fact, I 2 goes 
to zero as A —> 00 , so the full finite contribution is de¬ 
termined by I\. Therefore, the logarithmically divergent 
and finite terms coming from this part of IIv" (<z) are 


Hy,£;(g) 


9 2 
9 2 

kT 


27T 2 


27T 2 


87T 2 


ATn — 


IVln 


N In 2 


1 + z 2 
4 

1 + z 5 


In 



(61) 



1 

a + b cos 9 



cos 9 

a + b cos 9 


2 tt 

V a 2 - b 2 ’ 


27T 

T 


1 - 


\/ a 2 — b 2 


(53) 

(54) 


Doing these integrals, we find that Ily(g) splits into two 
parts, Ily(g) = II v,e{q) + n v,h(q)- These two parts are 


The remaining terms in the integrals in Eqs. (15311 and 
(l5il) give us YVv,h{q)- These are 

loi 2 a 2 r°° r 

f 00 , r , , sinh ii sin v sinh u! sin v' 

x / d/j, J dv - t - 

Jo Jo \ZF(p, is-,p',is') 











































X 


Ri — z 2 R 2 

( z 2 + cosh 2 g)(z 2 + cosh 2 y') ’ 


(62) 


0 


11 


where 

Ri = cosh y cosh y!{— sin 2 v sin 2 v' — \ cosh 2 y 

— \ cos 2 v — \ cosh 2 y' — i cos 2 v' + 1 

+ cosh y cos v cosh y! cos v), (63) 

i?2 = — cosh /x sin 2 v cosh y! sin 2 v' 

— ^ cos v cos v' (cosh 2 /i + cos 2 v + cosh 2 /i' 

+ cos 2 v' — 2 — 2 cosh y cos v cosh y cos v') , (64) 


F(/x, i/; y', i/) = [cosh(//+//) — cos(i/+ i/)] 
x [cosh(/x + y') — cos(i/ — i/')] 
x [cosh(/i — y') — cos(i/ + v')\ 
x [cosh(/x — y') — cos(f — i/')]. (65) 


This integral also has a logarithmic divergence. In fact, 
it too contributes to the In 2 term of IIy((?). Fortunately, 
we can determine this term analytically. We begin by 
introducing a pair of x and y coordinates corresponding 
to each pair of /i and v coordinates and by imposing a 
cutoff, as before. Doing this and simplifying the resulting 
expression, we obtain 


I q\ 2 o : 2 

n vM = ^^- 2 N I dy 

x / dy' I dx' ■ 


p2\+y 


J -1 

2A +y' 


dx 

1 


Rx - z 2 R 2 


\/ F {x,y,x',y ') 


(z 2 


>)(z 2 


- x' 2 ) ’ 


( 66 ) 


where A = 


v — Qo 
" «F|g|’ 


and F, i?i, and R 2 are F, R [, 
and R 2 , respectively, rewritten in terms of the a; and y 
coordinates. We now take the derivative of this integral 
with respect to A. The In 2 A term of the integral will cor¬ 
respond to a term in the derivative, with a coefficient 
twice as large as that of the In 2 term. The result that we 
obtain is 


3 , , 

d\ UvM 


I 

9 2 tt 2 

r 1 


f*2A -\-y 


= / d v 


dx 


dy' ■ 


1 


J-i \JF(x, y; 2A + y', y') 

R\{x' = 2A + y') — z 2 R 2 (x' = 2A + y') 


( z 2 + x 2 )[z 2 + (2A + y') 2 ] 


(67) 


We now expand the integrand in powers of j. The lead¬ 
ing term is, after integrating over y' , with respect to 
which the leading term is constant, 


d , \q\ 2 a 2 1 

d\ UvMq) " ^ 2 ^ N X, 

r 2X+v x 

/ dx —— 


dy 


■ 


( 68 ) 
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in 

II 

K? 

-20 
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-30 
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In A 


FIG. 5: Plot of for 2 = 5, where A = A/|<f| and 

z = qo/vF\q\- At large A, this curve becomes approximately 
linear; the slope of the line is exactly —4 and is twice the 
coefficient of the In 2 A term, while the constant term gives us 
the coefficient of the In A term. The red line is a linear ht to 
the points at large A. 


This is the only term in the expansion that will contribute 
to the In 2 A term of IIIt also contributes to the 
subleading In A term, and is in fact the only contribution 
with a coefficient that is a function of z; all other terms 
only contribute In A-divergent terms with constant coef¬ 
ficients. The total constant contributed by these, which 
we will call C , must be determined numerically. If we 
now evaluate the integrals, we obtain 


d\ 


n v,u(q) 


\qfa?_ 

g 2 2n 2 


N 



2 In A 


A 



(69) 


Therefore, the leading divergence of Uy,H (?) is 


UvMq) M ~^h N { ln2 {M 


+ 


In 


l + 2‘ 


c 


ln, il 


(70) 


The In 2 term exactly cancels that from the “easy” part, 
so that the whole diagram only has a simple logarithmic 
divergence. The full form of Uy (q) is thus 


n v (g) 


g 2 47r 2 


N 



—2(7111 



+ g(z) 


(71) 


where g(z) is a function that we will determine shortly. 

To obtain (7, we numerically evaluate A^-IIv ^(g) for a 
large value of A and z = 0 using Eq. dSZD and subtract off 
the In A term, which we know exactly. In the limit of large 


x 
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Z 

FIG. 6 : Plot of the function g(z) appearing in Eq. (1711) . giving 
the dependence of the finite part of IIv(<z) as a function of 
2 = qo/vF\q\- 


A, \ {q) is a linear function of In A, the linear term 

giving the coefficient of the In 2 term and the constant 
term giving the coefficient of the In term. We note that 
IIy ! _f/(g) and all two-loop contributions to the polariza- 

■ —m 2 2 

tion that we determine have the form, ^rj^NF(z,\), 
so that we only need to vary z and A. To perform these 
numerical calculations, we made use of the VEGAS algo¬ 
rithm as implemented in the CUBA package. We obtained 
results with errors of 0.1% or less; see Fig. [5] for an ex¬ 
ample, in which we plot A-^F(z,A) for z = 5. We find 
that 



FIG. 7: Two-loop rainbow correction to the electron self¬ 
energy, S 26 (?). 



FIG. 8 : Two-loop vertex correction to the electron self-energy, 
T, 2a (q). 


Straightforward algebra reveals that the integral over fco 
vanishes identically: 


dko 1 
27T k A 


7°(fco7° + VFk ■ 7)fc • 7(fco7° + VFk • 7 ) 7 ° 


= k • 7 


dk 0 kl~v 2 F \k \ 2 

2 ^(fc 0 2 + ^|£| 2 ) 2 


= 0 . 


(74) 


Therefore, the full contribution vanishes identically: 


S2 b(q) — 0. 


(75) 


C « 1.33318. 


(72) 


To obtain g(z), we evaluate Ih/,H(q) numerically for a 
large value of A and subtract off the logarithmically- 
divergent terms; our result is plotted in Fig. [ 6 ] 

IV. SECOND-ORDER SELF-ENERGY 


B. Two-loop vertex correction to self-energy 


Reduction to a quadruple integral. The second two- 
loop correction to the electron self-energy is shown in 
Fig. [ 8 ] This diagram has the value 


We now determine the electron self-energy to two 
loops. There are three possible diagrams that contribute 
to this order, and we discuss each in turn. 


A. Two-loop rainbow correction to self-energy 


£ 2a{q) 


-1 2 • 2 2 

— l07T ia v F 


d 4 k r d A p 1 1 

(2tt ) 4 7 (2?r ) 4 Iq’-kl 2 \q~p \ 2 


v _,o p _p r 1 r v r 0 

7 p 7 {k+p-qy 1 P 27 ' 


(76) 


The first potential two-loop correction to the electron 
self-energy is shown in Fig. 0 This diagram evaluates to 


The product of gamma matrices in the integrand can be 
expanded as 


^2 b(q) — 


(X 


/ d A k 1 

7°G 0 (fc)E 1 (fc)Go(fc)7° 

( 27r ) | q- k \ 2 


d A k 1 

( 2 n) A |g*_ fc |2 



7 ^' 7 ^ 7 


(73) 


yyyyyyy) = ^0 § u0 §P°sy 0 _ S v0 (1—S p0 )'y p 

—<5 M °(1—(5 iy0 )<5 p0 7 I/ —<5 M °(1—<5 1/0 )(1—5 p0 )7^7 p 7 ° 

— (1—(5 p0 )i5 iy0 (5 p0 7 p —(1—5 p0 )(5 1/0 (l— i5 p °)7 p 7 p 7° 

— (1—(5 p0 )(l—(5 l/0 )<5 p0 7 M 7 !/ 7° 

+(i-s il0 ){i-s v0 ){i-s p0 )^YY, 


(77) 
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which leads to 

S 2a (g)= 

= — 16TT 2 ia 2 v^ 


B 7 = 


' dk o f dp 0 p 0 


d 4 k f d 4 p 1 


' 2 n J 2 n M(k,p) 

1 Qo 


(2 tt ) 4 J (27r) 4 \q-k \ 2 W~p \ 2 k 2 p 2 {k+p-qY 

x ^k 0 p 0 (k 0 +p 0 -q 0 )'y Q -v F k 0 (k 0 +p 0 -q 0 )p-^ 

-v F k 0 Po{k+p-q}-^-v 2 F k 0 ^ 0 (k+p-q}-^p-^ 
-v F po(ko+po-qo)k-j-v F {ko+po-qoh°k-jp- ; Y 

—v%p 0 'Y°k-^(k+p—q)-'y+v%k-j(k+p—q)- ; yp- ;: Y (•. 


4v f \k\\k+p-q\[q 2 +v 2 F (\k\+\p\+\k+p-q\) 2 } 


(85) 


BfL = 


’ dfcp f dpo 1 


( 78 ) = 772 


27t J 2 n M(k,p) 

1 |fc|+|p|+|*+p-gl 


As usual, we first perform the integrals over the energies 
fcp and po : 


S, = 


S 2 = 


'dfco f dpp k 0 po(k 0 +po-qo) 

2 n J 2 n M(k,p) 

1_go_ 

4 go+ w F(l^l+|pl+l^+p-gl) 2 ’ 

f dko_ f dp 0 k 0 (k 0 +p 0 -q 0 ) 

J 2 n J 27r M(k,p) 
l |fc|+|p|+lfe+p-g| 

4 |p|[go+^F(l^l+|pl+l^+p-gl) 2 ] 


(79) 


(80) 


B = f dko_ f dp 0 k 0 p 0 


2 n J 2 n M{k,p) 

1_ \k\+\p\+\k+p-q\ _ 

4 \k+p-q\[qo+v 2 F {\k\+\p\+\k+p-q\) 2 } ’ 


(81) 


B = f dko f dp 0 k 0 


2 n J 2 n M{k,p) 

1 go 


4w f |p||fc+p-g|[gg+u|(|fc| + |p| + |A:+p-g|) 2 ] 


(82) 


S, = 


= 


’ dfcp r dp 0 Po(k 0 +po-q 0 ) 

2 TT J 2 7T M(k,p) 

1 |fc|+|p|+|fc+p-g1 

4 |fc|[go+ w F(l^l+|p|+l^+p-gl) 2 ] ’ 

’ dfcp f dp 0 k 0 +po-q 0 
2 TT J 2 TT M(k,p) 

_1_go_ 

Av f \k\\p\[q%+v 2 F (\k\+\p\+\k+p-q\) 2 ] 


(83) 


(84) 


Av f \k\\p\\k+p-q\[q^+v 2 F (\k\+\p]+\k+p-q\) 2 y 

( 86 ) 


with 


M(k,p) = 

(k^+v 2 F \k\ 2 )(pl+v 2 F \p\ 2 )[(k 0 +p 0 -q 0 ) 2 +v 2 F \k+p-q] 2 }. 

(87) 


We then have 


S 2Q (g)=-167r 2 ja 2 7;| J 


d 3 k f d 3 p 1 


(2 n) 3 J ( 2 n ) 3 \q-k \ 2 \q~p\ 2 
x |i?i7 0 — B 2 v F p-Ff—B 3 v F (k+p— g)-7 

—i?4 7 0 (fc+p—g) • 7p- 7 —S 5 uf • 7—7?6 7 0 fc • 7p -7 

-S 7 w|,7 0 fc-7(fc+p-g)-7+S 8 t;|,fc-7(fe+p-g)-7p-7|. 

( 88 ) 

The 7 0 component of this is 


i Tr[ 7 °E 2a ( (? )]= 

2 2 2 f d 3 k f d 3 p 1 1 r 

= -4 n 2 ia 2 v 2 F / ——r / — _ u q { 4-£?i 

J (2F) 3 i (2 f ) 3 |g*_jb|2 |g— p| 2 *■ 

—-4S4UFp-(fc+p— g)— ABeVpk-p—AB 7 v 2 F k-(k+p—q)} 

, 2-22 f d 3 k f d 3 p 1 1 

=~ 47 r 1 a v F qo 


1 


{ 2 tt) 3 J (2 tt ) 3 \q-k \ 2 |g—p| : 


1- 


q%+v 2 F {\k\+\p\+\k+p-q \) 2 
p k+p—q k p k k+p—q 


\p\ |fc+p—g| |fc| \p\ |£| \k+p-q\ J 


(89) 
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while the spatial components are 

1. 


-Tr[7 I E 2a (g)]= 


= —47 T 2 ia 2 Vp 


d 3 k f d 3 p 1 


-{-4 B 2 v FPi 


(2tt) 3 7 (2tt) 3 |g-fc|2 \q-p\ 2 
-AB 3 v f (ki +pi-qi)-AB 5 v f ki +4B$v F hp-(k+p-q) 

-( ki+pi-qi)k-p+pik-(k+p-q ) | 

d 3 p 1 1 


=—47r 2 ia 2 u^ 


' d 3 fc 

(2n) 3 J (2?r) 3 |£_fc|2 |<f-p| 2 


|fc|+|p| + |fc+p-g| / fcj Pi 

do+ v F (\k\+\p\+\k+p-q\) 2 \ \k\ \p\ 

h+Pi-qi k i (\p\ 2 ~p-q)+p i (\k\ 2 -k-q)+q l k-p 


\k+P~q\ 


\k\\p\\k+p-q\ 


switch to prolate spheroidal coordinates: 

\k\ 

\p\ = — (cosh p — cos v ), 

\k + p\ = (cosh p + cos v ), 

_ \k\ 2 

k-p = ——— (cosh p cos v — 1) , 

- \k\ 2 

p- (k+p) = —— (cosh 2 p + cos 2 v — 2) 

|fc| + |p| + \k+p\ = |fc| (cosh/x + 1), 

\k\ 3 


d 3 p = -sinh p sin ^(cosh 2 p — cos 2 u) . 

(93) 

Using spherical coordinates for k and performing the triv¬ 
ial integrations over 9 and over the angular variables as- 
(90) sociated with k. we find 


If we choose the coordinates such that <f=(0, 0, |g]), then 
it becomes apparent that the terms of the integrand 
which are proportional to k x ,k y or p x ,p y are odd func¬ 
tions of these variables, implying that these terms vanish 
upon integration. We may then make the replacement 
ki —> qik-q/\q\ 2 , and similarly for pi. The total integral 
is therefore proportional to qp 


-Tr[ 7 °E 2Q ( g )] = 


ia 2 v 2 F qo 


d|&||fc| / dpdv 


sinh p sin v 


irr, r w ( m -4tt 2 ia 2 v 3 F qi 
i Tr[ 7 S 2a («)]=-^- 

\k\+\p\+\k+p-q] 


+ 


d 3 k f d 3 p 


1 


(2tt ) 3 J (2tt) 3 \q-k\ 2 \q~P\ 
k-q p-q (k+p-q)-q 


sin 2 v sinh 2 % 

X ___ 

(cosh p — cosv) 2 ' 

Next we perform the integration over |/c|: 

1 


9o + v F \k\ 2 (coshp + l) 2 


(94) 


dl+v 2 F (\k\+\p\+\k+p-q\) 2 l |fc| \p\ \k+p-q\ 

\p\ 2 k-q+\k\ 2 p-q+\q\ 2 k-p—2(k-q)(p-q) 

|fc| \p\\k+p—q\ 


Tr[ 7 u E 2a (g)] = -4 tt ia v F q 0 
d 3 k f d 3 p 1 1 


(2t t) 3 J (2t r) 3 |fc|2 |^ 2 g 2 +u 2 (|fc| + |;p1 + |fc+^) 2 


1- 


p-(k+p) k-(k+p) k-p 1 
\p\\k+p\ \k\\k+p\ \k\\p\i 


(92) 


[ d\k\\k\ 
Jo 


do + v 2 F \k\ 2 (l + cosh/i) 2 


1 


2(1 + cosh p) 2 v F 


In 1 + 


do 


-(1 + cosh p)‘ 


(95) 


Off 


d define x = cosh p, y = cos v to obtain: 


Extracting the divergence in the temporal part:- To ex¬ 
tract the divergent logarithm term in the temporal part 
of the two-loop self-energy, Eq. (l89l) . we must exam¬ 
ine the behavior of the quadruple integral in the region 
|ft|, |p| » |g|. In this regime, the integral reduces to 


1tV[ 7 °E 2a (q)] = 


ia 2 go 

47T 2 

rl 


dx- 


x — 1 


dy 


'i ( 1+*) 2 

1 - y 2 


In 1 + 


(A 2 


?o 2 


(l + x) ; 


J -1 
ia 2 q 0 


(x - yy 


i: 


dx 


x — 1 


x (xarccothx — 1)] 


(1 + x) 2 


In 1 + 


?o 2 


-(l + x) ; 


ioi 2 qo (IQ _ t^ 2 ) ^ u f a 


27T 2 

ido 

2 n 2 


finite 


do 

A 


(10 — 7T 2 ) a 2 In — + finite. 
I <71 


(96) 


where we have discarded terms which become odd under 
the change of variable k —> —k, p —> —p in the limit 
of large |fc|, \p\. To perform the remaining integrals, we 


Extracting the divergence in the spatial part. To ex¬ 
tract the divergent logarithm term in the spatial part 
of the two-loop self-energy [Eq. (|91[) ]. we must exam¬ 
ine the behavior of the quadruple integral in the region 
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|fc|, \p\ |g|. It helps to first redefine k —> k+q: 

, v , \i -4? T 2 ia 2 v 3 F gi f d 3 k f d 3 p 1_f_ 

4 Tr [ 7 2a(g)] |gT JwJwWW 

y \k+q\+\p\+\k+p\ \ {k+q)-q p-q | (fc+p>g 

go+^l^+gl+lpl+l^+pl) 2 1 l*+g| IpI |*+p| 

| Ip1 2 (fc+g)-g+|fc+g1 2 p-g+|g1 2 fc-p-2(fc-g)(p-g)-|g1 2 p-g 

|*+g||p||*+p| 

(97) 

We then make the following expansions in the large mo¬ 
mentum limit: 


|/c+g|«|fc| + 


k-q 1 1 k-q 1 


1 p-q 


\k\ ’ \k+q\ \k\ |fc| 3 ’ \P~q\ |p| IpI 3 ’ 

\k+q\+\p\+\k+p\ 1 


q%+vl{\k+q]+\p\+\k+p \) 2 qo+v 2 F (\k\+\p\+\k+p \) 2 

go-^(l^l+|p| + l^+P|) 2 k-q 


\k\+\p\+\k+p\+ 


q$+ v 2 F (\k\+\p\+\k+p \) 2 \k\. 


= Q{k,p)+R{k,p) ^S, 
\k\ 


(98) 


where the last equality defines the functions Q{k,p) and 
R{k,p). We then have 


-2 — “ 


d 3 k f d 3 p Q{k,p) |p] ({k-q ) 2 |g| 2 


(2tt) 3 7 (2tt) 3 |jfc|2|p]2 |jfc+£| ^ |jfe|3 |jfc| ) ’ 

( 101 ) 


-3 — 


-4 = 


5=-2 


s =-2 


^=2 


d 3 k j 

d 3 p Q{k,p) {k-q){p-q) 

(2tt) 3 y 

( 2 tt ) 3 |fe|2|£|2 \k\\^\k+p\’ 

d 3 k I 

d 3 P Q{k,p) \q] 2 {k-p) 

(27t) 3 J 

(2^) 3 \k\ 2 \p \ 2 \k\\p\\k+p\' 

' d 3 k 

f d 3 p Q{k,p) ( k-q){p-q) 

(2t t) 3 J 

(2tt) 3 fc 2 \p \ 2 \k\\p\\k+p\ 

r d 3 k 

f d 3 p Q{k,p) (k-q){p-q) 

J (2tt) 3 

J (2tt) 3 |jfc|2|^|2 \k\\p \ 2 

o f d 3 k r d 3 p Q{k,p) {p-q ) 2 

“J (2tt ) 3 J (27 r) 3 |fc|2|^2 |pf ’ 

d 3 k r 

d 3 p Q{k,p) {p-q){k+p)-c, 

(2tt ) 3 J 

(2tt) 3 |jfc|2|^|2 \p\ 2 \k+p\ 


( 102 ) 

(103) 

=-2S 3 , 

(104) 

(105) 


(106) 


(107) 


iw[ 7 % 0 (g)] 


—4:TT 2 ia 2 Vpqi 

w 


d 3 k r 

WrJ 


X 



Q{k,p)+R{k,p) 


k-q 

w 


d 3 p 1 

( 2 tt ) 3 | jfe | 2|^|2 

p-q (*+£)•<? 

IpI \k+p\ 


{k+q)-q \ | \p\{k+q)-q L k-q \ 

1*1 v i*iv m+p\ v 1 * 1 v 

l*|p-g [ 1+ ^i ] | |g| 2 fc-p A *-q | 

IpP+pI v l*IV l*l|pll*+p| v l^lv 

2 {k-q){p-q) A fc-g \ |gl¥-q A *-g \ 1 

\k\\p\\k+p\ V \k\ 2 ) |*||p||*+p| V 1*1 2 /J 


(99) 


The terms that scale as the inverse sixth power in the mo¬ 
menta fc,p give rise to a logarithmic divergence. These 
are the terms we are interested in. There are also terms 
in Eq. HMD which scale as the inverse fifth power and 
so would seem to produce a linear divergence. However, 
these terms vanish identically as can be seen by perform¬ 
ing a coordinate transformation k —^ —k, p —> —p. We 
isolate each of the terms which contribute to the loga¬ 
rithmic divergence in the following series of integrals: 


d 3 k f d 3 p Q{k,p) {k-q){p-q) 


{2tt) 3 J (27t) 3 A: 2 |p| 2 \k p] \k+p\ 

o f d 3 k f d 3 p Q{k,p) \k\{p-q ) 2 

(2 tt) 3 j 

* (2tt) 3 fc 2 pj2 ^|3 ’ 

r d=k 

f d 3 p R{k,p) {k-q ) 2 

J { 2 n)- 

W (2tt) 3 |fc|2|^2 |fc|2 ’ 

f d 3 k 

r d 3 p R{k,p) (k-q){p-q) 

J { 2 n) 3 j 

' (2tt) 3 fc 2 p| 2 \k p| 

f d 3 k f 

d 3 p R{k,p) (k-q)(k+p)-q 

J { 2 n ) 3 J (27T) 3 |jfe|2|^|2 |fc||fc+p| 

f d 3 k 1 

r d 3 p R{k,p) \p\{k-q ) 2 

^ J (2 n) 3 J 

(2tt) 3 fc| 2 p| 2 \k 2 \k+p\ 


(108) 

(109) 

( 110 ) 

(111) 

( 112 ) 

(113) 


= f rf3fc f d3p Q$,p) f {k-q ) 2 \q \ 2 \ nnm 

11 J {2tt) 3 J (2^) 3 |fc|2| p1 2 ^ |fc|3 1 j 



d 3 k r 
(2tt) 3 y 


d 3 p R{k,p) (k-q) (p-q) 
(2^) 3 \k\ 2 \p \ 2 \p\\k+p\ ’ 


( 114 ) 
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To perform these integrals, we again make use of pro¬ 
late spheroidal coordinates for p and ordinary spherical 
coordinates for k , in which case we have 


Q(k,p) -> 

R(k,p) -> 


— |g|( r h '‘ +1) = 6(„, ifci>, 

gg+?;|,|fc| 2 ( co sh p+1) 2 

qo~v 2 F \k\ 2 {coshp+l) 2 _ ~ r 

- --- 72 = \ k \)- 

qQ+ y p\k\ 2 (cosh p+1) 2 

(115) 


In addition to the relations given in Eq. (1M1) . we also 
make use of the following: 



FIG. 9: Two-loop bubble correction to the electron self¬ 
energy, £ 2c (<?). 


^Tr[ 7 ‘£ 2a (q)] 


4--jriorvpqi (10 — 7r 2 )|g|^ 


\e 

«(10 — 7T 2 ) 
27T 2 


8 n 4 Vp 

2 1 ( A 
a v F qi In I 


In 


(~ 

Vl5l 

( 120 ) 


k-q = |/c||</| cos 6 k, 

1 * 1 , 


p ■ q = — [q x sinh p sin v cos 9 + q y sinh p sin v sin 9 


+ q z (cosh p cos v — 1 )]. 


(116) 


Here, 9k is the polar coordinate associated with k. After 
we plug these expressions into the S, ; , we first perform 
the integrations over 9 as well as over (f>k , the azimuthal 
coordinate associated with k. All the remaining four¬ 
dimensional integrals are then functions of q x and q y only 
in the combination q 2 + q 2 , which we may rewrite as 
\q \ 2 — q 2 . We then set q z = |^| cos 6 *^ and perform all 
the 9k integrations, which are easily done. Replacing 
cosh p —> x, cos v —> y, we then find the following results: 


10 


£* = 


15 


2=11 


dx dy d\k\ 


\e 

127T 4 j j 0 

\k\(l + x)(x + y - 2)(1 ~ y 2 ) 
[q 2 +v 2 \k\ 2 (l+ X ) 2 ](x-yr 
2 


£- = 




24?r 4 


dx J dy 


d\k\ 


\k\[q 2 -v 2 F \k\ 2 (l+x) 2 ](x-m-y 2 ) 


[q 2 + v 2 \k\ 2 (l + x) 2 ] 2 (x-y ) 2 


(117) 


Each 3D integral is logarithmically divergent, and the co¬ 
efficient of the divergence can be computed straightfor¬ 
wardly by first performing the |fc| integrals exactly and 
keeping only the term proportional to ln(t)frA/go)- The 
remaining integrals on x and y can also be done exactly. 
The results are 


10 

£* 

2=1 

15 

£* 

2=11 

and thus 


£ 


2=1 


(lO-TT 2 )!^ 2 /A\ 

12 ? t 4 v 2 f Uiy 

/m 

24t t 4 v 2 f V 1^1 y 


+ finite, 


+ finite, (118) 


fc^ln(£) 

8 n 4 v 2 F V l«| y 


+ finite, 


(119) 


Combining this with what we obtained for the temporal 
part, we then have for the full diagram 

S 2 a (g) = i 10 2n 2 (loT 0 + v F q- 7 ) a 2 In ^ ^ . ( 121 ) 


C. Two-loop bubble correction to self-energy 


The two-loop bubble correction to the self-energy is 
shown in Fig. [9] This diagram evaluates to 

S 2c (g) = -J ^£ 7 0 Go(<?-fc)7°nB(fc)H fc 2 , (122) 

where ^(< 7 ) is the result for the one-loop bubble dia¬ 
gram given in Eq. CD- To compute this integral, we will 
devise a streamlined approach to computing self-energy 
diagrams that contain vacuum polarization subdiagrams. 
In other words, consider the more general situation where 
ns(fc) is replaced by an arbitrary polarization function 
n(fc). Dimensional analysis tells us that n(fc) must have 
the following general form: 


n (k) 


J*l! 

V F 


F(z)lnl-4 I 

. vi*iy 


G(z) 


(123) 


where 2 = ko/(v F \k\) as usual, and F(z) and G(z ) are 
arbitrary functions. Using this form, we can rewrite our 
expression for the corresponding self-energy £(< 7 ) as 


A W> 


(124) 


We are only interested in the logarithmically divergent 
terms, so we expand the fermion propagator in the limit 
of large k: 


(q ~ k ) 2 


Ak? 


l-z 2 


(I + 2 2 ) 2 


<?o7 


z 2 + 1/3 


(1 + z 2 ) 


2 V VFq ‘ 7 


(125) 
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We then have 


V. RENORMALIZATION 


£(g) = - 


Mac 


dz 


1-z 2 
(1 + z 2 ) 2 


z 2 + 1/3 _ 

(TT^F ^ 7 


9o 7 


1 


/ ^\k\ ■ 
'151 \k\ 


+ G(z)\. 


F(z) In 



The integral over |fc| is easily performed, allowing us to 
express the self-energy as 

£(q) = | [A 2 g 0 7 ° - B 2 v F q ■ 7 ] In 2 

+ [Aiq-o 7 ° - Bxvpq ■ 7 ] In ^^ j , (127) 


where 


A 2 

B 2 

A! 

Bi 


dzF(z) 

-OO 

OO 

dzF(z) 

-OO 

dzG(z) 
dzG(z) 


1 -z 2 
(1 + z 2 ) 2 ' 
z 2 + 1/3 
(1 + z 2 ) 2 ' 
1 -z 2 


(1 + Z 2 ) 2 ’ 

Z 2 + 1/3 
(1 + ^ 2 ) 2 ' 


(128) 


We are now ready to evaluate £ 20 ( 9 )- From the ex¬ 
pression for ns(g) given in Eq. (Hill , we have 

N N /l + z 2 \ 

**> = -5?' = <129> 

Plugging these functions into the above expressions for 
Ai, etc. we find 

N N N 

A 2 -°, B 2 -- —, A!-- —, i? 1 -—. (130) 

These results then lead to the final expression for the 
self-energy correction: 

. . 21V 2 - _ 2 / A 

S 2c(9) = vpq-^ln 

+ 7^2 [ 3 ® 7 ° + VF $' fl “ 2 ln ( j|) ■ ( 131 ) 

Combining £ 2a and £ 2c , we then obtain the full 
second-order electron self-energy: 

v ( s 21V. 2 a /A 

S 2(9) = -^2«* v F q- 7hi 

15 + N 1\ o / 45 + IV 1 

+ (- 5 ^— 2 1 1 


x a 2 ln 


(132) 


We will now demonstrate the renormalizability of our 
theory and derive the renormalization group (RG) equa¬ 
tion for the interaction strength, Fermi velocity, and 
quasiparticle residue to second order in a. We will also 
derive identities that serve as valuable checks for our pre¬ 
vious results. 

All of the results derived in the previous section depend 
on an ultraviolet momentum cutoff scale A. In a complete 
theory (which is unknown) that takes into account the 
full band structure and interaction without any approxi¬ 
mation, the cutoff dependence would be absent since the 
energy dispersion will not be precisely linear all the way 
to large momenta in such an exact theory, and there¬ 
fore no explicit cutoff will be necessary. The cutoff indi¬ 
cates the energy/momentum regime where our effective 
low-energy theory defined by our action in Eq. ([ 6 ]) that 
describes the excitations at the Dirac cones is no longer 
valid. This high-momentum regime can be, for example, 
identified with the Brillouin zone boundary, the position 
of a van Hove singularity in the band structure, or simply 
the inverse lattice constant. We emphasize, however, that 
even such an unknown exact band theory, if it includes 
quantum interaction effects correctly, would manifest the 
ultraviolet renormalization effects in the sense that as the 
theory is used to calculate lower and lower momentum 
or energy properties of the system, it would manifest a 
logarithmic divergence as the Dirac point is approached, 
but the effective cutoff momentum entering such an ex¬ 
act theory would itself be also scale-dependent and not 
like the fixed cutoff A appearing in our formalism, which 
is based on a model linear Dirac band dispersion. In 
our effective theory framework, the cutoff separates the 
space of low-energy modes with known dispersion and in¬ 
teraction from the a priori unknown high-energy modes, 
and when restricting internal loop momenta to k < A, 
we essentially discard the latter. In a complete renor¬ 
malized theory, these contributions must be taken into 
account separately. For observables measured at small 
scales p -C A, the high-energy contributions arise from 
modes fluctuating over time and distance scales ~ 1/A, 
i.e., they appear as essentially local corrections. A sim¬ 
ple power-counting argument shows that only the self¬ 
energy, the vertex, and the polarization have a strong 
dependence on the high-energy cutoff, and we account 
for the missing high-energy contributions by adding local 
terms to the action. These terms are known as countert¬ 
erms. The coefficients of the counterterms are a priori 
unknown but should be chosen in such a way that our 
computation gives the experimentally observable renor¬ 
malized parameters, which only depend on the cutoff in 
subleading order 0(p/ A). As we discuss in the following, 
this procedure allows us to remove any strong cutoff de¬ 
pendence from the theory, yet introduces a dependence of 
the renormalized parameters on a low-energy renormal¬ 
ization scale. The divergence structure of the theory also 
determines the running of the renormalized parameters. 
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This section is structured as follows: First, in Sec. IV Al 
we introduce the renormalized Hamiltonian and discuss 
the relation between bare and renormalized parameters. 
This is followed by a detailed and abstract discussion of 
counter-term renormalization for 3D Dirac materials in 
Sec. IV Bl In particular, we derive various identities that 
link the coefficients of logarithmic divergences in self¬ 
energy, vertex, and polarization to the beta functions for 
coupling, Fermi velocity, and field strength. In Sec. IV Cl 
we use the results of the first part of this paper to derive 
the beta functions explicitly. The two-loop RG equations 
have some additional corrections stemming from the in¬ 
sertion of lower-order counterterms in the one-loop dia¬ 
grams, and we include a derivation of these diagrams in 
this section. Based on the results for the beta function, 
we comment on the critical interaction strength a c at 
which the RG flow has a fixed point, possibly, but not 
necessarily, indicating a breakdown of the perturbative 
expansion. 


A. Renormalized perturbation theory 

We first formulate a renormalized theory, in which we 
introduce counterterms that remove logarithmic diver¬ 
gences as a function of the momentum cutoff A from all 
results obtained from the theory. To do so, we relabel all 
fields and constants in the original theory with a super¬ 
script B to denote their “bare” values, and then intro¬ 
duce renormalized fields and constants. In other words, 
we rewrite the action as 

N 

s = -J2 

a—l 

+ ^b 7 °^b ) + _J_ J dt d 3 R( di (pB) 2 . (133) 

We first rescale the electron field, ip = Z^ X ' 2 ip B , and the 

Coulomb field, . If we perform this field 

rescaling and rearrange terms, the action becomes 

S B = S 0 + S cU (134) 


J dt d 3 R (^f + v F ip B Ydiipa 


call the renormalization scale. The second part, S c t, gives 
us the counterterms mentioned above. As discussed, they 
have a precise meaning in that they account for the high- 
energy physics that is not captured by the low-energy 
theory ©• The 5 coefficients in this part of the action 
are given by 


Sq 1 ^-0) 

(137) 

V B 

5i = l 

Vp 

(138) 

Sy = 1 - z^,zy 2 . 

(139) 

c _ , 9 1 y 

Op — 1 2 

(140) 


9 b 


The Feynman rules associated with So are the same as 
those given earlier for the “bare” action. Those for the 
counterterms are as follows. 


- X - 

- i(S 0 k 0 ^ a + v F Sik ■ 7) 

(141) 

a/x/XxX/x/v 


(142) 


0 r 

- Z7 S v 

(143) 


B. Identities 


As we have seen, when computing correlation functions 
defined by the theory given by Eq. (I133D . 

G %,m[{ x i}>{yjh<x B A] 

= (Tipa(x 1 )...ip B (x n )(p B (yi)...ip B (y m )}, (144) 

we will encounter divergences as a function of the mo¬ 
mentum cutoff A in our theory. The Green’s functions 
acquire an explicit cutoff dependence as indicated on the 
left-hand side. Our goal in defining the renormalized the¬ 
ory above is to eliminate this cutoff dependence, trading 
it for a dependence on the scale g. The renormalized 
correlators determined from this theory, 


where 

N 

So = -£ 

a—l 

+ ^a7 O, 0a) + ^2 J dt d 3 R (di(p) 2 , (135) 

N 

S« = £ 

a= 1 

+ Syipipa^fpa) - -^2 J dt d 3 R (diip) 2 . (136) 

The first part of the action, So, is just the “bare” action, 
but with a renormalized Fermi velocity v F and charge g 
that are defined at a momentum scale g, which we will 


Gn,m (aq,... ,x n \yi,... ,y m ;a,v F ,g) 

= Z~ n/ 2 G B ^x 1 ,...,x n \y 1 ,...,y rn \a B ,v%,A], 

(145) 

where 

Gn,m {{x t },{y t }]a,v F ,y) 

= (Tlp a (x l) . . . 1 p a (Xn)<p{y l) • ■ ■ <fi(y m )), (146) 

thus do not depend on the cutoff. The renormalized 
Green’s functions are instead functions of the renormal¬ 
ization scale g, the renormalized Fermi velocity v F , and 
the renormalized charge g 2 . As we will show in the fol¬ 
lowing, these renormalized quantities are running quan¬ 
tities that depend on the scale g, i.e., g 2 = g 2 (g) and 


dt d 3 R (lpa7°d t 1pa + Vplpa'fdilpa 


dt d 3 R (S 0 1pal 0 dt1pa + SiVpIpaYdilpa 
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vf = Uf (/-/)• with an evolution that is governed by renor¬ 
malization group equations. The cutoff dependence is re¬ 
moved by allowing the bare parameters (the bare Fermi 
velocity v B and the bare residue Z^f) to depend on the 
cutoff so as to cancel the divergences order by order in 
perturbation theory: 

Zijj = 1 - F^(a, A,p), (147) 

Zp = 1 — F v (a, A,p), (148) 

Vp = v F [1 - F vp (a, A, p)] , (149) 

(. g B ) 2 =g 2 [l-F 9 (a,A ,p)), (150) 

where 

Fj,(a, A,p) = ^ (a) In” (— \ , (151) 

n> 1 

F v (a, A, p) = ^ z v !n (a) In" f-'j , (152) 

F„ f (a, A, /x) = ^ u n (a) In” (■—\ , (153) 

n>l VA 4 / 

F s (a, A, /x) = ^ 5 „(a) In” . (154) 

The dependence on powers of In 0 f these functions 
is required to cancel the divergences in A coming from 
the bare action. The bare fine-structure constant a B can 
similarly be written as 

a B = a[l - F a (a, A,/x)], (155) 

where 

F a (a, A,p) = f n {a) In” f. (156) 

n>l vA 4 / 


The Callan-Symanzik equation then follows by using the 
chain rule: 


d d d 

'‘^ + / '"aS + ' ,F7 ~ay + " 7 * 

^ 5 i^F^p) = 0 , 

where we define 


da 


vf"1v f = ti¬ 


dy F 

dp 


7-0 y/z; — / 1 ^ \ Z -^0 
= A 4 ^ y/% 


(162) 

(163) 

(164) 

(165) 

(166) 


Note that, by construction , the n-particle Green’s func¬ 
tion G nim ({xj}; 5 ,uf,/ r) is finite, for we determine the 
coefficients in Eqs. (11471) . (11491) . and (11501) precisely to 
remove the divergent terms. This implies that f3 g , 7 „ F , 
and 7-0 must be free of any divergences in A. In addition, 
it turns out that /3 a , 7 „ F , 70 , and 7 ^ are independent of 
vf- 

Let us now relate p a , P VF , 7 ^, and 7 ^ to the coefficients 
of the divergent pieces in Eqs. (11471) and (I149D . First, let 
us consider the coupling a. Since the bare coupling a B 
cannot depend on p, 


p—a 

dp 


= 0 . 


If we now substitute in Eq. (11551) . we get 


(167) 


The functions f n can be expressed as a linear combina¬ 
tion of the various v n 's and e n ’s, i.e., 


/i(a) = 01 (a) - vi(a), (157) 

/ 2 (a) = 52 (a) - v 2 (a) +gi{a)vi(a) - fi(a) 2 , (158) 

/ 3 (a) = 53 (a) - v 3 (a) + 5 2 (a)ui(a) + 5 i(a)ui(a ) 2 
+ 5 i(a)u 2 (a) — 2 v\{ot)v 2 (pt) — Ui(a) 3 , 

(159) 


and so on. 

Let us now consider the bare n-particle correlator de¬ 
fined in Eq. (j!44l) . Note that by the definition of Sb in 
Eq. (11331) it does not depend on 5 , i.e., 

M^Gf,r”[{a ; J>{ 2 /i} ; a B ,u|, A] = 0. (160) 

We can use this fact to derive the Callan-Symanzik equa¬ 
tion for the n-particle correlator, Eq. (11461) . by combining 
Eq. (11601) with Eq. (I145D . This gives 




Gn^m ({^i}, a, Vf-j p) 


= 0 . (161) 


A^M 1 - F a (a,A,p)}} 


= 0 , 


(168) 


or 


Pa [1 - F a (a,A,p)] 

- pa4~[ 1 - F a (a, A, p)] 
dp 


= 0 . 


g B ,vp 


(169) 


We can rewrite this as 


Pa — F a (cq A, 5 ) Pa 

( d d d 

-a \ a ~da JrVFlvF ~fh>F + ^~dp ) Fa{°L, A, p). 


(170) 


Since P a must be finite, it is completely determined by 

/i(a): 


Pa=~ 


a/i(a). 


(171) 
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Higher-order terms obey the recursion relation, 


1. One-loop RG analysis 


(n + l)/„+i(a) = P a 


f n (a) + a 


df n 

da 


+ av F J VF 


dfn 

dvF 


(172) 


We can derive similar relations for the Fermi velocity 
anomalous dimension 7 „ F . Requiring that the bare Fermi 
velocity be constant in /.i yields 


lv F = 7 v f F vf (a, A, n) + 


d 


d 


VF ^d^ + ^ 


Pa-jr- 

F vp (a, A,/x), (173) 


from which, by a similar method as for g above, we obtain 
the relations, 


lv F = -vi (a) 


(174) 


and 


(n + l)v n+ i(a) = + Vf ^ Vf q^ 

x v n {a). 


(175) 


The divergent part of the first-order self-energy [Eq. 
©] is 


Si(g) = i^ln 



VFq ■ 7 + finite. 


(180) 


This divergence can be canceled out by the self energy¬ 
like counterterm (I141D . This, in fact, yields the values 
of 5 0 and to 0 (a ); we find that <5o = 0 (a 2 ), and 
thus there is no field strength renormalization to this 
order, and that 


^ = ^ ln (^)+ 0 (a 2) . (181) 

The leading-order polarization, as we see from Eq. CD- 
is divergent as well. This divergence may be canceled 
by the Coulomb propagatorlike counterterm, Eq. (114211 . 
This defines the renormalization S p to order O(a): 

= (182> 
Finally, the vertex correction at q = 0 is given by: 


We may also derive such a relation for the field anoma¬ 
lous dimensions 7 ^, (a) and 7 ^( 0 ). The relations for the 
renormalization of are 


7 */>(«) = 2 ^,i( a ) 


(176) 


and 


(wfl)^ 5 ji+l(®) — ^Afc(^) d^ ^7i/> Zijj,n (^) ■> 

(177) 


and those for ip are 


and 


Tv>{a) = 2 2 b ¥>(«) 


d 


(178) 


(n + l)z v , t nFi(a) = P a (a)~ - 27 ^( 0 ) ^,„(a). 


da 


(179) 


C. RG analysis 



0. (183) 


Any divergence here would be canceled by the vertex 
counterterm, Eq. (11431) . We see that there is no diver¬ 
gence, and thus S v = 0(a 2 ). 


2. Two-loop RG analysis 

The divergent parts of the two-loop contributions to 
polarization, given by Eqs. (El and CD , are 




2(z 2 + 4) 
z 2 + 1 


3 In 


z 2 + 1 


ln, l 


(184) 


and 


“ J f^ NCIn (i 


(185) 


We now pursue the RG analysis for our system, deriv¬ 
ing the counterterms and RG equations and performing 
some important checks on our results. 


The diagram corresponding to the counterterm that can¬ 
cels out the divergence in ns^(g) is shown in Fig. [TT3 and 
we will denote it by n SE,cr(q)- We find that it is given 
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FIG. 10: Diagram corresponding to the counterterm, 
II se,ct(<i), for the two-loop self-energy polarization diagram. 


by 



n sE,cr{q) 
= 2i5iN I 


d 4 k 

WT 4 


■tr 


• 0 1 ■ 0 1 7 -> l 

l l T/ - Tl v f^ ’ I7 


9 3tt 
\q \ 2 4a 2 


In A _ 


kl 


--In 

1 2 


2 o 2 Nhl ~ 
g 97r z /r 


ln il“ 


--In 

1 2 


z 2 + l 


(186) 


Since the vertex renormalization is zero to leading order 
(cf. Eq. (11831) 1. there is no contribution from 

= 0(a 3 ). (187) 


If we now add Ibs^ and Hse,ct, we find that it partially 
cancels the divergence from IIsb; all that is left behind 
are In 2 and In terms with constant coefficients: 



lls_B(g) + TlsE,cT(q) 

\e 2 a 2 


g 2 97t 2 


3 In 


N 


1 i 


2(z 2 + 4) 
z 2 + 1 


+ 2 + 1 
\e 2 cr 2 


9" 


97T 2 


N 


‘"'ll 

A 


- In ( - 
d 


- In ( - 

d 


(188) 


This remaining divergence plus that coming from Ily 
may then be canceled by adding the appropriate terms 
to S p . Upon doing so, S p becomes 




i + 


(V 8 

4^ 9 


a /A\ 2a 2 2 (^ 

7T J J 9 ?r2 \ M 

(189) 

We now turn our attention to the two-loop self-energy 
diagrams. There is only one nonvanishing second-order 
self-energy counterterm diagram, namely the one involv¬ 
ing the counterterm Coulomb propagator, shown in Fig. 
DU This diagram evaluates to 


L>2 ct — f<5: 


d 4 k 


+ 1 


( 2 tt ) 4 ' 7 (k 

47V o _ _ A A 

= d p T, 1 (q) = i—ra v F q • 7 In In —. (190) 
9 tt 2 |g| 9 


9? 


Y |fc| 4 g 2 


FIG. 11: (a) Coulomb counterterm contribution to the 

second-order self-energy counterterm diagrams, S 2 ct(q)- This 
is the only non-vanishing contribution, (b) Self-energy coun¬ 
terterm contribution. This diagram vanishes for the same 
reason that the “rainbow” diagram is zero, (c) Vertex coun¬ 
terterm contribution. This diagram only contributes at orders 
higher than a 2 . 


This has precisely the form needed to remove the In 2 
term from E 2 (q) provided we add appropriate second- 
order terms to 8 \\ 


c 2 a 

Sl = s 1 " 


A 


+ 


45 + AT 
97T 2 


2 N 
97r 2 


2 1 2 
i In 


1 


— 2 ' a ' ln 


A N 
/+ 

+ e>(a 3 ). (191) 


<5o also receives contributions at second order: 
'15 +AT 1 


<5o = 


37T 2 


2 ^ ln 


+ 0(a 3 ). (192) 


3. Beta functions and RG equations 


We are now in a position to determine the RG equa¬ 
tions describing the various couplings in our theory. We 
do this by determining the coefficients in the F functions 
defined by Eqs. (11471) (11501) . As we saw above, we can 
then immediately read off the beta functions and anoma¬ 
lous dimensions from these and thus obtain the associ¬ 
ated RG equations. 

Let us begin with Z^. From Eq. (11371) . we see that, to 
second order in a, it is just 

2* = l-(^- 4 )<* J l..(^)- (193) 

This immediately tells us z\^(a), 

~l+(ct) = y—^2 -2 ) a ’ ( 194 ) 


/15 + TV 

\ 67T 2 


and thus 7 ^,, 


= 


1 

4 


( 195 ) 
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The RG equation for Z ^ is thus 


and the latter is 


d In Zj/j 
din g 


15 +AT 1 
3n 2 


~ 2 a 


(196) 


Next, we will look at vj?. We first solve Eq. (11381) for 


»#. 


1 — <5i 1 — <5i 


vf 1 — < 5 o 

B 

Using our results, we find that, to order a 2 , ^ 

A 


(197) 


is 


u 


B 


2 a. 


= 1 -— In 
Up 3-7T 

2N 
97r 2 


2 i 2 

a In 


21V 

97T 2 


or In 


This immediately tells us v\(a) to second order: 

2a 2N 2 

n(a> = SF “ ■ 

The RG equation for up is then 

din up / \ 2a 2iV 2 

^ =7„ =-».(«) = ■ 


(198) 


(199) 


( 200 ) 


It is interesting to note that the second-order correction 
to the RG flow of the velocity only depends on the self¬ 
energy correction coming from the RPA-type bubble di¬ 
agram shown in Fig. [9l while the contribution coming 
from Fig. [ 8 ] cancels out exactly. 

Next, we look at Z v . As is shown in App. [GJ 1 — S v = 
Z^, and thus, to second order in a, Z tp = l. We therefore 
find that Z v remains constant for all values of /i. Finally, 

we consider g 2 . Solving Eq. (11401) for gives us 


{g B Y i 

g 2 1 - V 


( 201 ) 


Expanding to second order in a, we get 


(. g B ) 2 2 a 

~ = 1 + ^r N 

g 2 3?r 




2 a 2 

97T 2 


N{2N + 1) In 2 


In I — 


( 202 ) 


From this, we can now read off 51 (a): 


2 a 

9l(Q ) = --JV 


1 + l|c-?U 


(203) 


We can now determine /i(a) and thus the RG equation 
for a. The former is 


AO) 


-\(N+ 1)- 
3 7 r 



da 2(N + 1) 2 

= -aj i(a) =---a 


d In g 


37T 


27C - 44 . 
547T 2 


-Na d 


(205) 

These equations agree exactly with those that we would 
obtain directly from the Callan-Symanzik equation; we 
show this in Appendix [D] If we calculate the value of the 
coefficient of the a 3 term, then we find that it is negative: 


27C - 44 
54f 2 


-0.0150184 


(206) 


Therefore, we find that there are two fixed points for this 
equation, similarly to the case of graphene^, a = 0 and 


OLq 


36f 

44 - 27 C 



14.1298 



(207) 


If we start with a value below the latter critical value, 
a < a c , at a large scale g, then a will go to zero as we 
decrease g. On the other hand, if we start with a > a c , 
then we instead find that a diverges to infinity as g de¬ 
creases, indicating a runaway to infinite coupling. Similar 
behavior was also found in graphene^, although in that 
case the corresponding critical coupling is much smaller: 
a c = 0.78. In the present context of 3D Dirac materials, 
however, an additional interesting feature, not present in 
2D graphene, arises. Notice that the second-order correc¬ 
tion to the velocity RG equation, Eq. (12001) also differs in 
sign from the leading-order term, implying the existence 
of a second special value of a, this time at a* = 3n/N. 
If we start with a < a*, then the velocity grows mono- 
tonically as g decreases, while for a > a* the velocity 
is initially suppressed as g decreases. Moreover, since 
a* < a c for all values of N, there exists a window of 
couplings, a* < a < a c , in which both the coupling 
and velocity become smaller with decreasing g. This be¬ 
havior persists until a falls below a*, at which point the 
velocity reverses course and grows as the low-energy, non¬ 
interacting fixed point is approached. Thus, we find that 
for sufficiently large interaction strengths, the velocity 
exhibits non-monotonic behavior as the weak-coupling 
fixed point is approached at the Dirac point. We illus¬ 
trate this behavior in Figs. IT2l and fl3l in which we show 
a numerical solution of the RG equations for a and vf 
for a(go) = 0.9 and N = 12, typical values for the py- 
rochlore iridates. This nonmonotonic running of the ve¬ 
locity for a* < a < a c is in total contrast to graphene 
velocity renormalization, which is always monotonic for 
all coupling constants. 

One noteworthy feature of the RG flow for the Fermi 
velocity that creates a major contrast to graphene is the 
fact that, unlike in graphene, where the Fermi velocity 
diverges at low energ y 56 ! 60 : 65 (at least for a strictly non- 
relativistic, instantaneous interaction; in the real sys¬ 
tem with relativistic interactions, this growth will cease 
when the Fermi velocity reaches the speed of light) the 
Fermi velocity remains finite in the limit that the en¬ 
ergy scale goes to zero (i.e., at the Dirac point). This 
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FIG. 12: Numerical solution of the RG equations for a (red) 
and vf (blue) for a(fio) = 0.9 and N = 12, which are typical 
values for the pyrochlore iridates. Since a(fio) is far below 
a c , indicated by the red dashed line, we find that it increases 
for small /x, then saturates to a c . vf, on the other hand, 
decreases as /x increases up to about p, « 0.95/xo, then starts 
increasing. This is better illustrated in Fig. 1131 below. 



0 


FIG. 13: Plot of vf alone, with the range for p adjusted to 
better illustrate the non-monotonic behavior of vf- 


can be traced back to the fact that, in graphene, charge 
does not renormalize, as shown in Table [I] As a result, 
since graphene becomes weakly interacting at low-energy 
scales, the Fermi velocity must diverge. This is not the 
case for 3D Dirac materials. Because charge can renor¬ 
malize in this case, there is another way for the system 
to become weakly interacting at low energy, and that is 
for the charge to go to zero. This is exactly what hap¬ 
pens here. As a result, the Fermi velocity remains finite, 
even at low energies; it simply renormalizes to a different 
value. Thus, the nature of ultraviolet renormalization in 
3D and 2D Dirac materials is qualitatively different. 

As in the case of graphene^, the appearance of a c and 
ck* in 3D Dirac systems may be a manifestation of the 


breakdown of perturbation theory (due to the asymptotic 
nature of the perturbative expansion in a) rather than 
real quantum phase transitions, particularly since these 
points precisely correspond to the second-order correc¬ 
tions becoming comparable in magnitude to the leading- 
order results. Following Dyson’s original argument for 
QED2£, we can obtain an estimate for the order of the 
expansion beyond which the perturbative results can no 
longer be truste d 65 ’ 84 : 


n 


3.09797 

a 3 / 2 


(208) 


Here, n is the expected perturbative order up to which 
the loop expansion in the effective fine structure con¬ 
stant is asymptotic, with the perturbation series di¬ 
verging beyond the nth order. As an example, let us 
consider CdsAs 2 : the experimental values^ - — of the 
Fermi velocity are typically in the range 10 5 — 10 6 m/s, 
while the dielectric constant has been measured to be 
k ~ 20 — 40 ) 67 i 68 suggesting coupling values in the range 
a ~ 0.1 — 1. Equation (12081) then suggests that our 
second-order perturbative analysis is likely to be valid for 
coupling values toward the lower end of this range, but 
as a approaches unity, the results become questionable. 
If we interpret a* as a signature of the breakdown of the 
perturbative expansion, then we would further conclude 
that the validity of perturbation theory also depends on 
N since a* is inversely proportional to N. For a mate¬ 
rial like CdaAs 2 for which N = 1, we find a* = 37r ss 9.4, 
which is more or less consistent with Eq. (12081) . However, 
for the pyrochlore iridates^ with N = 12, we instead ob¬ 
tain a* = 7 t/ 4 « 0.8, which may call into question the 
applicability of perturbation theory for these materials 
if their effective interaction strength a satisfies a > ck*. 
The fact that perturbation theory would seem to be more 
reliable in the regime of smaller N is interesting given 
that this is precisely complementary to the regime in 
which RPA or a large -N expansion is vali d 62 ’ 66 . It would 
appear that perturbation theory and RPA/large-iV are 
in this sense complementary to one another. 

In contrast to QED, where the bare coupling is fixed 
to be ~ so that the perturbative expansion is asymp¬ 
totic to a very high order, the 3D Dirac materials enable, 
in principle, an investigation of both the running of the 
coupling for a specific material and the analysis of the 
asymptotic nature of the underlying field theory by vary¬ 
ing the Dirac material so as to modify the bare coupling 
from a value as small as 0.1 to as large as 1. 


4- Consistency check 

We will now apply the recursion relations for the coun¬ 
terterms, Eqs. (11721) and (11751) . to our results to ensure 
that they are satisfied, thus providing a compelling check 
on our calculations. To be exact, we can show that our 
results for f\(a) and V\ (a), at order a, along with the 
beta functions and anomalous dimensions determined in 
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the previous section, give the same values of f 2 (a) and 
^ 2 ( 0 ) to order or that we determined directly from our 
two-loop calculations. 

Let us start with the relation for a, Eq. (I172D . We 
already found fi{a) in the previous section, so we only 
need to determine / 2 (a). Using the previously-stated 
formula for this function, we find that, to order a 2 , 


f 2 (a) = g 2 {a) - v 2 (a) + gi(a)vi(a) - vx{a ) 2 

= -^( N + l) 2 . (209) 


Now let us check that we obtain this result from our 
recursion relation. If we substitute /1 (a) and into Eq. 
(11721) and only keep the lowest-order terms, we obtain 


2 0/2(0:) 


2(iV + 1) 2 
3 tt “ 


1 + a 


da 


8 (N + l) 2 3 
97r 2 


2(N + 1) 
37r 


a 


( 210 ) 


This gives us the value of f 2 {a ) stated above. 

We may similarly apply the relation for vf , Eq. (11751) . 
We may immediately read off v 2 {a) from our previous 
results: 


■*(“) = 9 ^ ( 211 ) 

Again, we wish to show that we obtain this same re¬ 
sult from the recursion relation. Substituting Vi (a), 7 „ F , 
and /3 a into Eq. (11751) and keeping only the lowest-order 
terms, we find that 


2 v 2 {a) 


(2a 2(N + 1) 2 d \ 2a 

\ 37T 37t da) 37t 

4AT , 

—-pra . 


( 212 ) 


We thus obtain the same value of v 2 {a) stated above. 


VI. EXPERIMENTAL DETECTION 


While initial experiments have reported control over the 
bulk detuning by varying the surface doping^, this might 
not yet be sufficient to provide definite signatures of 
renormalization. In addition, typical samples can have 
strong disorder. However, a clear observation of scaling 
does not require the suppression of all relevant perturba¬ 
tions, but rather a hierarchy of scales where one parame¬ 
ter (for example, the temperature) can be tuned to large 
values compared to any other scale (for example, the dop¬ 
ing). In this section, we discuss this largely overlooked 
alternative way of detecting renormalization effects by 
varying the temperature. 

This section is structured as follows: First, we give 
an example for the integrated charge renormalization to 
leading order in a which can be obtained in closed ana¬ 
lytical form, thereby introducing the Landau pole scale 
which signifies a divergence in the renormalized coupling. 
Second, we illustrate the effects of finite temperature on 
the Drude weight as obtained from kinetic theory as a 
generic example of an observable that can be detected 
in optical conductivity measurements. We show that the 
Drude weight of an extrinsic system (i.e., with finite dop¬ 
ing) at finite temperature assumes the same value as the 
Drude weight of an intrinsic system, i.e., increasing the 
temperature probes the intrinsic system irrespective of 
initial detuning. The renormalization is then apparent in 
an additional temperature dependence of the quasiparti¬ 
cle parameters giving rise to a deviation from the intrin¬ 
sic unrenormalized linear-in-temperature scaling. In the 
third subsection, we present a calculation of the plasmon 
behavior to leading order in RPA, which naturally incor¬ 
porates the effect of finite doping and the renormalization 
of a. Thus, by studying plasmon properties experimen¬ 
tally as a function of doping density and/or temperature, 
it should be possible to verify the ultraviolet renormal¬ 
ization and the RG flow of 3D Dirac system many-body 
effects. 


A. Coupling constant renormalization and Landau 

pole 


In this section, we discuss how the running of the quasi¬ 
particle properties discussed in the previous sections can 
be detected in experiment. Such an experiment has to 
extract either the Fermi velocity or the fine-structure con¬ 
stant over a parameter range to check for scaling. Previ¬ 
ous experiments on 2D Dirac materials had precise con¬ 
trol over temperature, doping, and disorder, all of which 
are competing scales that cut off the intrinsic renormal¬ 
ization group flow at low energies, possibly masking the 
intrinsic Dirac point physics. For example, varying the 
doping by means of an external gate potential, the Fermi 
velocity renormalization induces an increase in the Fermi 
velocity as smaller densities (corresponding to smaller 
Fermi energies) are probed^. 

It appears that at the present time, a comparable con¬ 
trol over 3D Dirac material properties is still lacking. 


To illustrate the coupling constant renormalization, we 
consider the renormalization group flow of the coupling 
to one loop. As derived in Eq. (12051) . the renormalization 
group equation is 

da _ 2(N+l)a 2 
d In g 37t 

If the beta function has negative sign, the coupling be¬ 
comes weaker with increasing energy scale - this is the sit¬ 
uation encountered in quantum chromodynamics, where 
the perturbative high-energy limit is known as ‘asymp¬ 
totic freedom’. In the present theory, however, the cou¬ 
pling increases with increasing energy scale, just like in 
quantum electrodynamics. Starting in a perturbative 
regime at small energy, the system becomes nonperturba- 
tive at higher energy. Solving the renormalization group 
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equation (12131) with the boundary condition a(fio) = ao 
(the measured value of the charge at a scale fi 0 ), we ob¬ 
tain: 




OtQ 


2(N + l)o!o fi o 

1 +---In — 

67T fl 


(214) 


It turns out that this solution has a divergence at a scale 
fj, = Ap = /j J0 e 3 ' K / 2< ' N+1 ' >a ° known as the Landau pole. 
We can express the renormalized coupling (12141) in terms 
of the Landau pole: 


a(n) 


3tt ... i A L 
2{N+1) n fi 


(215) 


The Landau pole is thus the effective parameter that 
characterizes the interacting Dirac semimetal. At scales 
close to Ax,, the coupling diverges and the theory be¬ 
comes strongly interacting, even if we start with a weakly 
interacting theory at low energy. The theory makes no 
prediction for Ax; it is a material property that must be 
taken from experiment. Strictly speaking, as the renor¬ 
malized coupling diverges, our first-order theory loses va¬ 
lidity, and a more complete and fully nonperturbative 
calculation is required in order to establish if the diver¬ 
gence is cut off by nonperturbative effects or if the cou¬ 
pling truly diverges. Our theory, however, makes the 
unambiguous prediction that the charge should increase 
logarithmically at high energy or temperature (setting 
fi = T). By analogy with quantum electrodynamics, one 
would expect that the pole is preempted by a phase tran¬ 
sition, although this remains to be established for Dirac 
semimetals^. 

We mention that by definition our leading-order beta 
function [Eq. (I213D ] has no asymptotic freedom (i.e., 
there is a positive sign on the right hand side for all a val¬ 
ues); in fact, we have the Landau pole at very high energy 
where the running coupling diverges. We note, however 
(as discussed in Sec. [Q of our paper), that our calculated 
beta function to second order [c.f. Eq. (12051) ] does have 
a negative sign if a > a c with the renormalized coupling 
running to infinity at lower-energy scales (i.e., behaving 
superficially similar to the QCD situation) provided we 
start in the strong-coupling regime of a > a c . In such 
a strong-coupling situation (which arises only when the 
second-order perturbative correction is comparable to the 
first-order term, and hence may very well be unreliable 
within a loop expansion theory), the Dirac system man¬ 
ifests stronger (weaker) coupling at low (high) energy, 
superficially mimicking QCD behavior. 


B. Optical conductivity and Drude weight 

In this section, we show the competing effects of tem¬ 
perature and doping on the Drude weight of a 3D Dirac 
material and demonstrate that even a moderate temper¬ 
ature of order T ~ £p/ 2, where ex is the doping-induced 
Fermi energy, is sufficient to probe the intrinsic response 


cs 
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FIG. 14: Drude weight for extrinsic (red line) and intrinsic 
(green dashed line) Dirac materials. At high temperature, 
there is no difference between the intrinsic and the extrin¬ 
sic Drude weight, thus manifestly showing that temperature 
can act as an appropriate running scale to study the intrinsic 
Dirac point physics. 


of the system without detuning. The Drude-Boltzmann 
form of the low-frequency conductivity takes the form 


cr(w) 



47 T UJ — i/r' 


(216) 


where w 2 is the Drude weight and 1/r is a transport re¬ 
laxation time induced by disorder or interactions. Within 
kinetic theory, the Drude weight is given by (a detailed 
derivation is given in Appendix lEl) 


w 2 d = 4?re 2 


2v 2 p Df 

“4T 



de J’ 


(217) 


where Dp = N £ 2 F / [k 2 v'p) is the density of states and /o 
the Fermi-Dirac distribution. The red curve in Fig. 1141 
shows the extrinsic Drude weight as a function of tem¬ 
perature. The low-temperature Drude weight (which is 
dominated by intraband excitations) starts off at a con¬ 
stant value with a Sommerfeld correction that decreases 
its value at small but finite temperature. The Drude 
weight assumes a minimum and then increases at high 
temperature where it has a linear temperature depen¬ 
dence. For comparison, we include a plot of the intrinsic 
Drude weight as given by Eq. (12171) with the chemical 
potential set to zero. It is apparent from the figure that 
at high temperature both curves are equal and the intrin¬ 
sic Drude weight is measured, even if the system started 
off at finite doping. The renormalization should induce a 
superlinear high-temperature dependence reflecting the 
intrinsic semimetallic many-body effects as discussed in 
the next section. 

The computation of the Drude weight within kinetic 
theory gives the same results as a diagrammatic calcu¬ 
lation using the Kubo formula to the leading one-loop 
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order—. Note that in the intrinsic limit, higher-order 
diagrams are expected to introduce nonlocal corrections 
to the conductivit y 87 ’ 88 . Here, we consider the extrinsic 
case, for which this effect does not appear. 


C. Plasmon dispersion 


We now consider the plasmon dispersion which is given 
by the pole of the dielectric function. It is related to the 
Drude-Boltzmann form of the conductivity by 


e(w) = 1 — = 1—^=0. (218) 

A full solution of the plasmon dispersion to leading order 
in RPA was presented in Ref. [H, and will not be repeated 
here. Instead, we expand on this work by considering the 
effect of charge renormalization. The defining expression 
for the dielectric function does not involve the bare cou¬ 
pling a, Fermi velocity, and cutoff separately, but only 
combined in the form of the Landau pole without any ad¬ 
ditional divergent terms. The full solution of Eq. (12181) 
then gives the plasmon dispersion as a function of tem¬ 
perature, Fermi energy, and Landau pole. The plasmon 
dispersion is shown in Fig. [15] for both doped extrinsic 
(red line) and undoped intrinsic (green line) Dirac sys¬ 
tems using a Landau pole scale of Al =8 ep. As dis¬ 
cussed in the previous section, doping does not influence 
the plasmon at high temperature and the dispersion is 
equal for intrinsic and extrinsic configurations. The full 
calculation of the plasmon dispersion which includes the 
coupling constant renormalization predicts a superlinear 
high-temperature scaling of 




T 

In A l /T 


a(T)T. 


(219) 


This superlinear scaling is valid at high temperature and 
is cut off a low temperature due to finite doping, as seen 
in Fig. 1151 While simple dimensional analysis would only 
predict a linear temperature dependence of the intrin¬ 
sic plasmon,— the correct result displays a superlinear 
scaling through the logarithmic temperature dependence 
due to the renormalization of the coupling as discussed 
in Sec. IVI Al Note that the result for the intrinsic plas¬ 
mon agrees with the Drude weight (12171) if we assume 
a constant k and replace the bare value of charge and 
Fermi velocity by the corresponding renormalized expres¬ 
sion (12151) . 

An experimental observation of such a super-linearity 
in the measured temperature dependence of the plasmon 
dispersion would be strong evidence for the ultraviolet 
renormalization of interaction effects in Dirac materi¬ 
als. We emphasize that since 3D Dirac materials can¬ 
not be doped in situ in a continuous manner using an 
external gate (as can be done for 2D graphene layers), 
temperature-dependent measurements (in the regime of 
temperature being larger than the Fermi energy associ¬ 
ated with the doping density) are probably the simplest 



FIG. 15: Plasmon dispersion for extrinsic (red line) and in¬ 
trinsic (green dashed line) systems with A l = 8 ef- The 
manifest superlinearity in the temperature-dependent plas¬ 
mon dispersion of the doped system arises from the ultravio¬ 
let renormalization with temperature also acting as the cutoff 
for T -< Tf, as expected. 


technique to characterize their intrinsic Dirac point prop¬ 
erties. Such temperature dependence will have a low- 
temperature cutoff (or saturation) due to the remnant 
doping in the sample, as shown in Fig. [151 but above 
this cutoff scale of the extrinsic Fermi energy, the system 
should reflect the RG scaling defined by its Dirac point 
properties. 


VII. CONCLUSION 


We have calculated the electron self-energy and polar¬ 
izability for Dirac-Weyl materials up to second order in 
the effective fine-structure constant a and for arbitrary 
number of pairs of Dirac cones N (arising, for example, 
from valley degeneracy associated with band-structure ef¬ 
fects). As expected, we find that log-squared divergences 
arise at second order as necessary for the renormalizabil- 
ity of the theory, and we have used these to perform 
many self-consistency checks of our results. After cal¬ 
culating the counterterms necessary to cure these diver¬ 
gences, we extracted from them the many-body renor¬ 
malization effects that are evident in properties such as 
the quasiparticle residue, Fermi velocity, and coupling 
strength. Our derivation of the RG flows of the quasi¬ 
particle properties reveal the existence of a critical point, 
a c = 14.13(1 + 1/N), above which the coupling flows to 
larger values at low energies. This behavior is very analo¬ 
gous to what occurs in graphene, which also exhibits such 
a fixed point, albeit at a much lower value of the coupling: 
a c = 0.78^. Following Dyson’s original argumen t 65 ' 83 , 
we have argued that the appearance of a c may be an ar¬ 
tifact of the asymptotic nature of the perturbative series 
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and not a real feature of the system, but only measure¬ 
ments can determine whether there is a strong-coupling 
fixed point in 3D Dirac systems. It is conceivable that our 
perturbative expansion is correctly indicating the pres¬ 
ence of a strong-coupling fixed point for the system, but 
the precise value of the critical coupling strength can only 
be correctly evaluated by non-perturbative methods. 

Unlike graphene, however, Dirac-Weyl materials ex¬ 
hibit a second special value of the coupling, a* = 3n/N, 
at which the flow of the velocity reverses direction. This 
implies that for couplings in the range a* < a < a c , 
both the interaction strength and Fermi velocity initially 
decrease with the energy scale until the coupling reaches 
the value a,, at which point the velocity flow reverses 
and begins to climb to larger values as the energy scale 
is reduced further. At this point, it is unclear whether 
a* should be viewed as giving rise to an experimentally 
observable non-monotonic behavior in the Fermi velocity, 
or if it should instead be seen as an indication that per¬ 
turbation theory cannot be trusted for coupling strengths 
comparable to or larger than a* (as is likely to be the case 
for a c ). This is a very interesting avenue for further ex¬ 
ploration given that a* may be in the experimentally rel¬ 
evant regime for materials possessing a large Dirac cone 
multiplicity, such as the pyrochlore iridates. We hope 
that future experimental work in 3D Dirac materials will 
look for possible signatures of nonmonotonicity in the 
measured Fermi velocity as a function of doping density 
(at low temperatures) or temperature (at low doping den¬ 
sities) to check for the predicted existence of a* in the 
theory. 

We also discuss the experimental detection of ultra¬ 
violet renormalization effects by computing the Drude 
weight and the plasmon dispersion up to leading order 
in the fine-structure constant and the random phase ap¬ 
proximation, respectively. There are two main findings: 
first, by performing computations at finite doping and 
temperature, we establish that the intrinsic renormaliza¬ 
tion group results which we obtain in this work govern 
the behavior of the system at high temperature, i.e., the 
behavior of a Dirac material is the same as for an intrin¬ 
sic system irrespective of initial detuning. This general 
result (that holds not only for the observables discussed 
in this paper) is of particular relevance for current ex¬ 
periments which are not able to precisely control doping 
and disorder. Second, we find that the high-temperature 
Drude weight including renormalization and the plasmon 
dispersion scale in a superlinear form as ~ T / In A j j /T. 
While the term linear in temperature can be predicted 
from dimensional analysis, the logarithmic scaling is a 
direct manifestation of renormalization. The strength of 
the renormalization effect is set by the Landau pole scale 
Al, a scale which resembles similar findings in QED at 
one-loop order— and which provides a renormalization- 
flow invariant parameter that may be used to fit our re¬ 
sults to experimental data. The superlinear temperature 
scaling should provide a direct signature of renormaliza¬ 
tion effects in current and future experiments on Dirac 


materials. 
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Appendix A: First-order polarization bubble 


Here, we give the details of how to determine the first- 
order polarization bubble n B (g), shown in Fig. [T] The 
expression that we obtain is 


n B (q) = N j ^Tr[ 7 0 Go(/c + g) 7 °Go(fc)] 


= -N 


drk 

( 2^) 4 


Tr 


(k + q ) 2 ^ k 2 


.(Al) 


Evaluating the trace, we obtain 


n B (g) = -AN 


d A k 


ko(k 0 + g 0 ) - vjpk ■ (k + q) 


(2ti-) 4 (fcg + v 2 F \k\ 2 )[{k 0 + g 0 ) 2 + v 2 F \k + g| 2 ] 

(A2) 

Evaluating the integral over fco, we get 


n B (g) = -2N J 


d 3 k 


\k\ + \k + q\ 


1 - 


^ ( 2 j) 3 q 2 Q +v 2 F (\k\ + \k + q\) 2 
k- (k + q) 


|fc||fc + g| 


(A3) 


If we now switch to prolate spheroidal coordinates, as 
given in Eq. (12211 , we obtain, after performing the (trivial) 
6 integration, 


n B (g) = - 


W-n[ 

8 tt 2 vf Jo 


cosh p 


z 2 + cosn p 


dp / dv sinh p sin v 
Jo 

(1 — cos 2 i/). 


(A4) 


We now make another coordinate change, x = cosh p 
and y = cosza Note, however, that the integral has a 
logarithmic divergence, and thus we introduce a cutoff 
|k| < Aon the momentum. Doing these, we obtain 


n B (g) = - 


kT 


(*2A +y 


87 t 2 vf 


N 


dy 


dx ■ 


r(l -V 2 )- 
(A5) 


If we now evaluate the x integral, we get 

n B (g) = - N J dy {l- y 2 ) In - 


+ (2A + y) 2 
z 2 + 1 


(A6) 

If we now perform an integration by parts, we obtain Eq. 
m in the main text. 
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Appendix B: First-order electron self-energy 


The diagram corresponding to the first-order electron 
self-energy correction is shown in Fig. [3j and will be de¬ 
noted by T,i(q). The expression that we obtain for this 
diagram is 

r h 

Sr(g)=- J ^ 7 °Go(Hg) 7 °flo(-fc), (Bl) 


or, upon substituting in the expressions for the bare 
Green’s function Go(fc + q) and the bare Coulomb prop¬ 
agator, 


£i(g) = -*5 2 


d 4 k 0 ft + ftl 0 1 
(2tt) 4 7 (.k + q ) 27 |fcja' 


(B2) 


The integral on fco may easily be done; after evaluating 
it, we obtain 


£ 1 ( 9 ) = 


d 3 k 0 {k + q)-j 0 
-IT -r7 ■ 


( 2^) 3 ' \k\ 2 \k + q\ 

Since 7 0 anticommutes with 7 *, we may rewrite this 


(B3) 


as 


„ , A I .2 f d 3 k (fc + oO-7 , , 

Ei(g) = 7*9 / 77-77 ^—-■ (B4) 

2 J (2tt) 3 Ifcplfc -)- gl 


We now note that the integrand is symmetric under ro¬ 
tations of k about the axis along q. Because of this, we 
may make the replacement, 


,T , , {k + q)-q^ „ 

(k + q) ■ 7 -A-r -77 - q ■ 7 , 


thus obtaining 


£i(g) = i 


kT 


g 2 f d 3 k (k + q) ■ q 


2|?1 2 J (2tt) 3 |fc| 2 |fc + g| 


9'7- 


(B5) 


(B 6 ) 


We now switch to spherical coordinates, with the +z axis 
in the direction of q. The integral now becomes 


£ 1 ( 9 ) = i 


2 I 9 T ( 2 tt )2 y 0 


d\k\ / sin# 


I 9 I 2 + |fc||g| cos$ 

|fc | 2 + \<f\ 2 + 2 |fc||q| cos 0 


=9 • 7- (B7) 



FIG. 16: Two-loop self-energy vertex diagram, Vse(p, q)- 


Appendix C: Second-order vertex diagrams 

In this appendix, we compute the two-loop corrections 
to the vertex and show that Z v = 1 identically to second 
order. 


1. Two-loop self-energy correction to vertex 
function 


The first diagram that we compute is the two-loop self¬ 
energy diagram, Fig. 1161 It is given by 


Vse(p,9) = 


d 4 k 


( 27 t ) 4 "* 7 ft-f ' ft 


*7°-^£i(fc)-p*7° 


\k - p + 9| 2 

(Cl) 

We isolate the logarithmic divergence by setting p = 0: 




Expanding the gamma matrix part, we find 
2 g 2 avF . f d 4 k 1 


(C2) 


1 ; S . e ( 0 , 9 ) ~ -- 


(3v F k%\k\ 2 - vf.|fc| 4 ) 7 0 


3tt J (2 tt) 4 k 6 \k + q\ 2 


(* 


+ (fco - 3u B fc 0 |fc| 2 ) k • 7 


4 1 A 

7 +ln — 
3 1*1 


(C3) 


It turns out that the fco integration vanishes identically, 
and thus VsE{p,q) is ultraviolet finite. 


Here, we impose a cutoff A on the |fc| integral because 
the integral is divergent. Naively, we would expect this 
integral to have a linear divergence, but we will see that 
it is only logarithmic. Let us make the substitutions, 
|fc| = \q\n and x = cos 9. The integral now becomes 


£1(9) 




1 + KX _ _ 

/, - =9-7- 

V 1 + K 2 + 2nx 


(B 8 ) 

These integrals can be performed exactly; doing so, we 
obtain Eq. in the main text. 


2. Two-loop vacuum polarization bubble correction 
to vertex function 


Next, we consider the vertex function diagram with 
the vacuum polarization bubble inserted in the internal 
Coulomb line, Fig. [T71 and we again set p = 0 to extract 
the divergence: 


Vb(0,<j) = 


d 4 k . 0 i . 0 i . 0 9 

11*7 7*7 7*7 


( 2 tt ) 4 ' ft 


\k + q\' 


-n B (fc + q). 


(C4) 
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FIG. 17: Two-loop vacuum polarization bubble vertex dia¬ 
gram, V B (p,q). 


FIG. 19: 
Vcc(q)- 



Two-loop crossed Coulomb line vertex 


diagram, 



FIG. 18: Two-loop parallel Coulomb line vertex diagram, 
Vpc{q )■ 


Since the divergent part of does not depend on fre¬ 
quency, the integration over ko for this term will be iden¬ 
tical to what we had for the one-loop vertex function 
diagram [Eq. (11831) ] . and it will thus vanish. Thus, the 
only potential divergence comes from the finite term in 

n B : 


is 


Vpc(p,q) = 


d 4 k 


d 4 £ 


( 2 tt ) 4 J ( 2 tt ) 4 * 7 tf + f $ 

. n i 


7 l l 


xi 7 


7l 7 


o g 2 g 2 

\k\ 2 \£\ 2 ' 


+ / 

(C 8 ) 


We isolate the potential logarithmic divergence by setting 
p = q = 0 in the integrand. If we then redefine k 
k — £, then it becomes straightforward to perform the 
integration over ko'. 


'dk 0 ft£t = / = 0 


fc 4 


k 4 


(C9) 

We conclude that this diagram has no ultraviolet diver¬ 
gence. 


4. Two-loop crossed Coulomb line correction to 
vertex function 


VsM 


ig 4 N 0 f d 4 k k% — v F \k\ 2 ^ 1 + z 2 


12 tt 2 v f J (27t) 4 k 4 \k\ 2 


(C5) 

Rewriting ko in terms of z and performing the trivial 
integration over fc, we find 


V B (0 ,q) 


r ~- ^- 1 .lnl±i!ln A 


- 7 U / dz- 


48tt 5 Vp 4 q\' 


(C 6 ) 


The integral over z evaluates to 7 r, leaving us with 

V B (0, q) = i-^a 2 7 ° In + finite. (C7) 
3 tt 2 |g| 


3. Two-loop parallel Coulomb line correction to 
vertex function 


We now consider the diagram with two parallel 
Coulomb lines, Fig. [TS] The expression for this diagram 


Here, we consider the diagram with two crossed 
Coulomb lines, Fig. [T9] The expression that we obtain is 


Vcc(p,q) = 


x vy 


d 4 k r d 4 


ij° ——7*7° 


(27t) 4 J (2t r) 4 4 + t d + ^ + 


p+f+f+t fl + j + ft \k\ 2 \£\ 2 


i ■ o 9 2 g 2 




(CIO) 


We again set p = q = 0, redefine £—>■—£ and then 
k —> k + £ to obtain 

d 4 k f d 4 i ^ 0 jt-f°h°h 0 (/ + W 




(27r) 4 pk 4 (k + £) 2 \k + l\ 2 \£\ 2 ' 


(Cll) 

Expanding out the gamma matrix product and throwing 
away terms which integrate to zero, we are left with 

„ , , ■ 4 o /' d 4 k f d 4 £ 

VccM~-*g 7 Jj^J W) 4 

x-—-—— ((kn - v 2 F \k\ 2 ){£l - v 2 f \£\ 2 ) 

£ 2 k 4 {k + £) 2 \k + £\ 2 \£\ 2 >■ 
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FIG. 20: Two-loop vertex correction vertex diagram, Vv{q)- 


- 3k 0 £ 0 v%\k\ 2 +klt 0 - 


VcoM~-w£rf$ s f 


(| k\ + \£\ + \k + £\r 


1 + 


k-i 

w 


( 2 tt ) 3 \I\>\k + i\* 

I • (C13) 


Next, we switch to prolate spherical coordinates in k , 
yielding 


Vcc(p,q) ~ -n 


o 5 


d 3 £ 1 


167t 2 ?;|, 


dpdv 


(27t) 3 |£|3 „ 

. . 1 — COS ^ 

x sinh p sin v- - -—— 

r (1 + coshp) 2 

cosh /i cos v + cosh p — cos v — 1 
(coshp + cos v) 2 

Performing the integrations over £ and v , we find 


Vcc q) ~ i7 ° 1 / ^ (l+cthri 2 


(C14) 


x (coshp — 1)(1 + coshp In tanh ^). 


The p-integral evaluates to — 


10 — 7T 


giving 


0 10 — 7T 2 , 

Vcc(p, <?) = -*7 9?r2 « In j-7 + finite. 


(CIS) 


(C16) 


5. Two-loop vertex correction to vertex function 

Finally, we consider the two-loop vertex correction di¬ 
agram, Fig. [20] This diagram gives us 


C d 4 k r d 4 £ 

VvM = 1 dk I de 


(27T) 4 J (270^ 0 + ^'' {! + # + ; 

i . r 




X Z 7 


ji 7 


2 2 

* . o 9 y 

j'i'y 


+ ft ]j) -\- (/[ + f |fc| 2 |£| 2 


^7° 


d 4 fc /• d 4 f tl°(f + lt)7 0 h 0 h 0 


(2tt) 4 7 (27t) 4 £2( fc+ ^2 fc 4|£|2|^2 ■ 

(C17) 


We proceed by observing that 

7 °(^ + = 7 °h°h°(tt + th° + ■■■, (cis) 

where the ... represent terms which vanish upon inte¬ 
gration over k and £. This means that the product of 
gamma matrices evaluates to precisely what he had in 
the previous subsection: 


4fco4> + 3/cg — v F \k\ 2 Vpk ■ /j . 

(C12) 

Performing the integrations over ko and £o, this becomes 

[ _d?k_ f d 3 £ 1 

Up 1 f ^ 3 1 

\k\ + \£\-\k + £\ 


V v {p,q) ~ *5 4 7°y 

x |(fcg - Vp\k\ 2 )(il - Vp\£\ 2 ) - 3k 0 l 0 v 2 p\k\ 2 


(2tt) 4 J (2tt) 4 pk 4 {k + £) 2 \k \ 2 \£\ 2 


+ kll 0 - |4fc 0 4 + 3 k 2 - 4|fc| 2 J v 2 pk ■ f} , (C19) 


and the integrations over fco and £q are also identical: 

• 4 r d 3 k f d 3 £ 1 


W(M)~n°!j -J 


U F J (27T) 

|fc| + |ij-|fc+ii 


3 I (2tt) 3 |fc| 2 |^| 2 |fc + £| 

k-e\ 


i + 

(|fc| + |^| + |fc + £|) 2 v 

Again using prolate spherical coordinates, we find 


(C20) 


Vv(p,q) ~ *7° 


5 


d 3 ^ 1 


167t 2 i;|, 


(2tt) 3 |4 3 „ 

1 — cos ^ cosh p cos 1 / + cosh p — cos — 1 


dpdv sinh p sin;/ 


(1 + coshp) 2 


(coshp — cosi') 2 


(C21) 


This is nearly identical to the analogous expression for 
Vccip, q)i except for the overall sign and the sign in the 
denominator of the final factor in the integrand. The 
latter comes from the fact that now we had |fc| 2 in the 
denominator instead of \k + £\ 2 as before. Despite this 
sign difference in the denominator, the integral over v 
evaluates to precisely the same value as before, so that 
the only difference compared with VcciPi d) is the overall 
sign: 


, s . o 10 - tt 2 2 , A 

Mp,4) = >7 -f-s-o, ln^ 


finite. 


(C22) 


It would then seem that the divergences of Vcc(jP, d) and 
V v (p, q) would cancel each other. However, Vy(p,d) re¬ 
ceives an extra symmetry factor of 2. 


6. Two-loop vertex function 

Summing up the various contributions to the vertex 
function, we find the total result 

v (2) (p, q) = V B {p, q ) + Vcc{p , d) + 2 V v (p, q) + finite 
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N 10-7t 2 \ , A , 

h, m + (C23> 

This implies that the relevant counterterm coefficient is 

^ = ( 3 ^ + “ 2 hl D + ° {a3) =So + ° {a3) - 

(C24) 

As pointed out in the main text, this implies that Z^ = 
1 to second order in a. This confirms the requirement 
S v = S 0 stemming from gauge invariance. 



Appendix D: Callan-Symanzik RG analysis 

Here, we present an alternative means of deriving the 
RG equations for our system, namely through direct ap¬ 
plication of the Callan-Symanzik equation. We will find 
that our results agree with those derived in the main text, 
providing a valuable consistency check. 


1. Photon two-point function 


The Callan-Symanzik equation for the photon two- 
point function expanded to third order in a reads 


+ (3 a d a T VF^fypdyp T 27 ^) 

v q a w 2 n(°) + v„ 2 nw + w 3 ( n (o)) 2 l = 0 . (di) 


Here, V q = ^TravF/\q\ 2 is the Coulomb interaction, and 
n(°) and n« are the leading- and next-to-leading-order 
vacuum polarization functions. Summing up the relevant 
diagrams and their respective counterterms from Secs. Mil 
and [V] yields the following expressions for the renormal¬ 
ized vacuum polarization functions: 


(0) = _A^ / A \ 

6^v f Vki; 


j% 1 2 

127 T 2 Vf 


In 



n (1) = 


a 


Na\q\ 2 
IStt 3 ^ 
2z 2 1 
1 + z 2 



Na\q\ 2 
18 tt 3 vf 
Na\q\ 2 
87 t 3 vf 



kl 

(D3) 


In tf 1 ), the term involving the constant C comes from 
the two-loop vertex correction to the polarizability. All 
other terms in come from the self-energy correction 
to the polarizability. 

We can use Eq. ED to derive equations for j3 a , 7 VF , 
and 7 ^ by plugging in n*°' and n^ 1 ) and solving the 
resulting equation order by order in a. The first-order 
equation we obtain is trivial: 


CD 4tt av F _ 

Ha |^j 2 — U ’ 


(D4) 


and simply tells us that the first-order beta function for 
a is zero: /3^ = 0 . 

The second-order equation reads as 

V q 2 fid^+^d a V q +v F 'yi 1 Jd VF V q +2 1 ^V q = 0. (D5) 

Plugging in the expressions for V q and n (0 \ we obtain 
the following relation: 

-«7$ ~ 2a 'r!p ) - ( D6 ) 

Next, we consider the third-order equation which fol¬ 
lows from Eq. ED: 

+ k, 3 ^(n (0) ) 2 + pWd a v q 
+/ 3 i 2 ) a ct (v g 2 n (0) ) + vp^JdypVg + VFi$d VF {v 2 n^) 
+2^ 2 W q = 0. (D7) 

Plugging in n(°\ n^ 1 ), f3a \ 7 ^, we find that the z de¬ 
pendence cancels out, and we are left with the following 
relation: 

Pa ] + + 2 a7^ 2) - [c - = 0. (D8) 

The fact that the z dependence cancels out of the Callan- 
Symanzik equation constitutes a nontrivial consistency 
check of our results for the vacuum polarizability. 

In order to obtain explicit expressions for /3„ ’ and , 
we must first calculate the first- and second-order contri¬ 
butions to 7 „ f and 7 ^ from the Callan-Symanzik equa¬ 
tions for the electron two-point function and electron- 
photon three-point function. 


2. Electron two-point function 


We can obtain 7 Vp as well as the scaling function for 
the electron field strength, 7 ^, from the Callan-Symanzik 
equation for the electron two-point function expanded to 
second order: 


(/r<9/i -|- (3 a d a A vf'Yvf d VF A 2y^/?) 
= 0 , 


-Si-Ei- 


-A-Et-A-22- 
i i $ i i 

(D9) 


This equation involves the renormalized first- and second- 
order electron self-energies, which we obtain by summing 
up the results from Secs. m and m 

El = ^ ivp w 7 ln ■ (D10) 


v 21V . 2 -z I 

L 2 = —-—^ria Ui?(7'7ln — 


97 T 2 


kl 
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15 + iV 1 

3tt 2 2 

45 + iV 1 

9tt 2 2 


907 


v F q ■ 7 


a 2 In 



(Dll) 

Plugging these results into the first-order equation that 
follows from Eq. (ID91) . 


3. Electron-photon three-point function 

To complete the analysis, we must also obtain 7 V , the 
scaling function for the photon field strength. We can 
compute this from the Callan-Symanzik equation for the 
electron-photon three-point function. Summing up the 
various contributions to the second-order vertex function, 
we find the renormalized result 


(vf + 7^) l ~ + = 0, (D12) 

and using the identity, 

= (dib) 

t i i 

yields the following leading-order behavior of j VF and 7 ^,: 

7j )= 0. (D14) 

We now turn our attention to the effects of the two- 
loop corrections to the electron self-energy. Collecting all 
the second-order terms in the Callan-Symanzik equation 
gives 


V< 2 >( M ) = * 7 ° (£5 + ^+) « 2 >" f + 

(D20) 

while we found in Sec. [V] that the first-order vertex func¬ 
tion has no divergence, V^{p,q) = finite, and is thus 
independent of the renormalization scale. The second- 
order Callan-Symanzik equation reads as 


(P^d a + v Fl <gd VF + 7 ™ + 27 ( ^) r(p, q) 




V p U^(p)T(p,q) + V p l -V^\p,q) 


+T(p,q)Ei(p + q)- 


Si(g)r(p, q) 


= 0, 
(D21) 


-ud^ 1 - + 


+^' > ^d a ^ 1 1 - + + VFjipdvp 


W-. 


+2^ = 0. 


(D15) 


Using the identity from Eq. (ID13I) and collecting all the 
finite terms, we find 


15 + iV 1 \ i 


r 2 


+9o7 


45 + IV 1 

9t r 2 2 


x -iv F q • 7- + 7^ ~iv F q ■ 7- + 27^ - = 0. D16) 

I $ $ i w $ 

Writing 


6 0 

- = —-zg 0 7 “ - ~i v fQ * 7-7 > 


(D17) 


and collecting like terms, we obtain the following two 
relations: 


a 


37T 2 
45 + iV 1 

9t r 2 2 


15 + iV 1 
2 


- 27?> = 0, 


+ 7$-27?’ = 0, (D18) 


which imply 


( 2 ) 2 /15 + iV 1 

7V = « '- 


67T 2 


4/ ’ 


~( 2 ) = 
IVf 


2Na 2 
97r 2 


(D19) 


where V), = 47rauF/|p| 2 , and where we have defined 


r (p,q) = V p l -i7° 

i 



(D22) 


Equation (ID21I) can be drastically simplified by using 
the Callan-Symanzik equations quoted earlier, namely, 
Eqs. (1D12I) and (1D5I) . yielding 


(-7^ 1 ) -2 7 J ) )r(p,g) = 0, (D23) 

which implies 

7^ = -27? } = 0. (D24) 

Inputting this result and the result for 7 ^ from 
Eq. (1D14|) into Eq. (1D6I) then produces an explicit ex¬ 
pression for the leading-order beta function for a: 


fl(2) = 2 ( 7V + 1 ) q 2 

Pa 37t a 


(D25) 


( 2 ) 

To obtain 7 ^, ; , we need to consider the Callan- 
Symanzik equation at third order in a. The full form 
of this equation is very complicated; however, we may 
make use of Eqs. (lD7l) and (ID 151) to simplify it consider¬ 
ably, with the result 

(-7^ 2) - 2 q. 0 2) ) *7° + ( /?i 2) da + v F 7vJd VF ) V ( 1 ) (p, q) 
+pd p V {2 \p,q) = 0. (D26) 

Although we have not directly computed V ' 1 - 1 (p, q), we 
can see from this equation that it must be independent 
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of momentum since there is no other momentum depen¬ 
dence in the equation. This would mean that (p, q) = 
V*V(0,0). However, we have shown in Eq. (11831) that the 
latter vanishes, hence, V^(p, q) = 0. Using Eq. (ID20I) . 
we then have 


7 


( 2 ) = _ 
V 





10 — 7T 2 \ 
2tt 2 ) 


a 2 = 0 . 


(D27) 


(2) (2) 

Since we now have 7 ^ and 7 ^ , we can go back and 
determine pa' 1 from Eq. (ID 8 I) : 


A 3) = -ajg ~ 2a 7 J, 2 > + (c - | 

(44 - 27 C)Na 3 
547r 2 


(D28) 


We thus find that we obtain all the same results as we 
did in Sec. [V] of the main text. 


Appendix E: Drude-Boltzmann form of the 
conductivity 


In this appendix, we derive an expression for the Drude 
part of the conductivity without renormalization effects 
using kinetic theory. Denote the distribution function 
momentum q at a node a with chirality s by f a s{k : t). 
In the absence of an external field, the noninteract¬ 
ing distribution is given by the Fermi-Dirac distribution 
fasio) = [exp[(siip’|g| — /i)/fc B T] + l] _1 . The full distribu¬ 
tion in the presence of an external electric field exerting 
a force F = h= eE on the electrons solves the Boltz- 

at 

mann equation 


l + P-E 

8 t dq 


fas(q,t) = - 


fas{q,t) - f° s (q,t) 


(El) 


where we assumed a relaxation time ansatz for the colli¬ 
sion integral. In linear response, the distribution function 
can be expanded as 


fas(q,u) = 2t r5(uj)f° s (q) + g as (q,uj), (E2) 


where g(q, cj) is determined from the Boltzmann equa¬ 
tion: 


VjJ - 


Q J 


ffas (g,w) = ^fas(q) 


Hence, 

ffas(g,w) = 


1 “I 


_eE d 0 
h ' dq fas ^ 


(E3) 


(E4) 


and the current is (assuming a constant scattering time) 


The conductivity is 


cr(w) = - 


1 


2 2 

v F e g 


ho+1/t 3tt 2 (Hvf) 3 Jo 
This yields the dielectric function 


dee 2 I — 


d f°(e) 
de 


UJZ, 


471-7 

g(w) = -cr w) = 2 ■ > 

i0 — UdJ 


where 7 = 1/r and w 2 denotes the Drude weight 


(E 6 ) 


(E7) 


2 2 d f f c 

= ^Tre —— —2~ / 

J e F do 


dee 2 - 


d f°(e) 
de 


(E 8 ) 


Here, Dp = N e 2 f / [k 2 (hv f) 3 ] is the density of states. 
Aside from the temperature dependence of the energy 
average, there is an additional temperature dependence 
of the chemical potential, which starts at g, = ef at zero 
temperature and then falls off to zero at high temper¬ 
ature. The extrinsic chemical potential g(T) at finite 
temperature is fixed by the relation 



d 3 k 

( 2 ^ 


n+(fc) + 2 N 


d 3 k 

(2^)3 


[n-{k) - 1 ] 


n, (E9) 


where n s are the Fermi-Dirac distributions with energy 
e s (k) = shvpk and n = gk F /6n 2 is the zero temperature 
carrier density. Equation (IE9I) reduces to 



(E10) 


The asymptotic behavior of g is 


ef 




1 



2 


T <C T F , 


T»T+. 


(Ell) 


This behavior is illustrated in Fig. [51] 


The temperature dependence of the Drude weight (|E 8 I) 
is determined both by the energy average and the tem¬ 
perature dependence of the chemical potential. At small 
and high temperature, we obtain 


dM = v F e^2 


d 3 q q 

w 1 ! 


fas(q,Ul) 


1 v F e 2 g 


iuj + 1 /r 3tt 2 (Hvf) 3 Jo 


dee 2 


de J 
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FIG. 21: Chemical potential vs. temperature as obtained by 
solving Eq. (IE10I) . Dashed lines indicate the low- and high- 
temperature limits (|E11I) . 

In the Born approximation, the transport scattering time 
is given by 

/q(fc) ~ f°a(k) 

T k 

= -^p / ^[mv\k')\ 2 fa{m~fa{k')) 

- \{k'\V\k)\ 2 f a (k')(l - /„(£)) }<% - e p ) 

= ~9a(k) j ^3^-ejf/) |(fc|^|fc / )| 2 (l-cos0) 

(E13) 


and hence 



(E14) 


Note that this transport scattering time is different from 
the quasiparticle lifetime derived from the imaginary part 
of the self-energy to this order in the interaction by the 
additional factor— of (1 — cos 9). In the present case, we 
assumed that Tg is energy-independent. 
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